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CHAPTER  ONE 


INTRODUCTION 

Progress  in  the  prediction  of  electromagnetic  interference  (EMI) 

* 

has  occurred  on  many  fronts  in  recent  years.  A  particularly  important 
development  in  EMI  prediction  involves  the  nonlinear  modeling  of  the 
Bipolar  Junction  Transistor  (BJT)  and  computer-aided  analysis  proced¬ 
ures.  A  large  number  of  pioneering  research  efforts  have  cul¬ 

minated  in  the  successful  EMI  analysis  of  the  bipolar  integrated 
circuit  (IC)  operational  amplifier  (Op  Amp).^^  The  engineering  ef¬ 
fort. and  computer  time,  and  hence  the  cost,  required  for  EMI  analysis 
of  the  full  IC  Op  Amp  are  large.  The  high  costs  involved  in  predict¬ 
ing  EMI  using  a  full  IC  Op  Amp  model  may  affect  the  wide  application 
of  the  powerful  computer-aided  EMI  analysis  techniques  which  have  been 

r 

developed.  The  research  which  will  be  described  in  this  report  uses  a 
simplified  model  called  a  macromodel  for  EMI  analysis  of  Op  Amps.  In 
using  the  Op  Amp  macromodel  for  EMI  analysis,  the  situation  of  the 
practicing  electromagnetic  compatibility  (EMC)  engineer  is  considered. 
An  Important  goal  of  the  present  research  is  that  techniques  general¬ 
ly  not  available  to  the  EMC  engineer  (such  as  the  laboratory  character 
ization  of  Integrated  electronic  devices)  are  to  be  circumvented. 


This  report  Is  organized  in  the  following  manner.  Chapter  One 
contains  a  review  of  material  relevant  to  this  report.  An  introduction 
to  the  theory  and  techniques  of  analysis  of  a  weakly  nonlinear  system 
is  presented.  Next,  an  analysis  of  EMI  in  the  741  IC  Op  Amp  made  with 
a  full  model  is  reviewed. Finally,  the  IC  Op  Amp  macromodel  de¬ 
veloped  by  Boyle  et  al .  is  presented  and  discussed. In  Chapter 
Two  the  nonlinear  macromodel  for  the  BJT  IC  Op  Amp  derived  from  the 
Boyle  model  is  presented.  Predictions  made  using  the  Nonlinear  Circuit 
Analysis  Program  (NCAP)^9^  and  the  nonl inear. macromodel  for  the  741 
Op  Amp  are  compared  with  full  model  predictions  and  experimental  results. 
Substrate  effects  on  macromodeling  accuracy  are  also  presented.  In 
Chapter  Three,  nonlinear  macromodeling  techniques  are  presented  in  detail. 
The  LMIO^57^  IC  Op  Amp  is  used  as  a  concrete  example.  It  is  shown  how 
the  macromodel  BJT  parameters  and  circuit  parameters  are  extracted  from 
published  information.  The  NCAP  predictions  on  how  amplitude  modulated 
RF  signals  are  demodulated  in' an  LM10  buffer  amplifier  are  compared  with 
experimental  results.  In  Chapter  Four  the  network  theory  for  the  non¬ 
linear  macromodel  is  developed.  Based  on  the  nonlinear  transfer  function 
approach,  it  is  shown  that  an  IC  Op  Amp  which  consists  of  cascaded  stages 
with  internal  feedback  has  nonlinear  characteristics  dominated  by  the  in¬ 
put  stage.  The  analysis  shows  that  the  input  BJTs  which  are  in  two  cascade 
configurations  can  be  represented  with  a  pair  of  single  BJTs.  Chapter 
Five  is  the  conclusion.  The  effectiveness  and  validity  of  the  non- 
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linear  macromodel  is  discussed.  A  nimber  of  new  research  topics  arising 
from  this  research  are  suggested. 

1.1  EMI,  Nonlinearity,  and  Nonlinear  Transfer  Function. 

As  communication  channels  become  increasingly  crowded,  unwanted 
signals  often  accompany  the  desired  signal.  The  unwanted  signals  can 
cause  interference  which  can  be  classified  as  being  caused  by  either 
linear  effects  or  nonlinear  effects.  Linear  Interference  effects  may 
occur  when  the  frequency  of  the  interfering  signal  is  very  close  to 
that  of  the  desired  signal.  Nonlinear  interference  effects  may  occur 
when  the  unwanted  signals  cause  new  signals  to  be  generated  by  the  non- 
linearities  in  the  electronic  devices  and  the  new  signal  frequencies 
are  close  to  that  of  tne  desired  signal.  Linear  interference  effects 
can  often  be  avoided  by  improved  tuning  and  selectivity.  The  nonlinear 
interference  effects  are  more  complex  and  not  subject  to  any  universal  ’ 
solution. 

Among  the  nonlinear  electromagnetic  interference  effects,  inter- 
modulation  is  an  important  example.  Intermodulation  interference  may 
occur  when  signals  with  frequencies  outside  of  a  receiver's  pass-band 
are  combined  nonlinearly  to  produce  signals  with  new  frequencies 
which  then  fall  within  the  receiver's  passband. 

Intermodulation  EMI  is  a  result  of  electronic  device  nonlinearity. 
As  illustrated  in  Figure  1-1,  two  sinusoidal  signals  with  frequencies 
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f1  and  f2  falling  outside  the  passband  of  a  receiver  having  second  de¬ 
gree  nonlinearities  produce  signals  with  new  frequencies. some  of 
these  new  frequencies  may  fall  within  the  receiver's  passband  and  be¬ 
come  interfering  signals.  The  nonlinearities  being  referred  to  in 
this  work  are  termed  "weak"  as  opposed  to  the  switching  type  of  non¬ 
linearity  found,  for  example,  in  a  modulator. 

In  order  to  analyse  an  electronic  circuit  to  predict  EMI,  equiva¬ 
lent  circuits  are  used  in  which  the  electronic  devices  are  represented 
by  nonlinear  models.  The  BJT  can  be  represented  by  the  Nonlinear  In¬ 
cremental  T-model  shown  in  Figure  1.2,  for  a  device  exhibiting  weak 
nonlinearities  under  small-signal  operation. £17*22^  in  this  model,  the 
nonlinear  base-emitter  resistor  is  represented  by  *he  nonlinear  incre¬ 
mental  current  generator  K(v2).  The  nonlinear  collector  capacitance 
Is  represented  by  the  nonlinear  incremental  current  generator  yc(v3-v2). 
The  dependent  current  source  in  the  collector  circuit  containing  the 
hp£  nonlinearity  and  the  avalanche  nonlinearity  is  represented  by  the 
nonlinear  current  generator  g(v2,v3-v2).  The  nonlinear  emitter  cap¬ 
acitance  is  represented  by  the  nonlinear  incremental  current  generator 
Ye(v2).  The  resistors  rfa  and  rc  are  the  incremental  base  and  collector 
resistances  respectively.  For  the  case  of  low-level  interference  (usual 
ly  labeled  as  weak  interference),  the  nonlinear  incremental  current 
generators  are  represented  by  Taylor  series. This  is  a  basic  dif¬ 
ference  between  the  small -signal  linear^-45^  and  nonlinear  BJT  models. 


i  I 


g(v,.v,-v,) 


y.  v,-v,) 


figure  1-2.  A  Small-Signal  Nonlinear  Incremental 
T-model  for  the  BJT. 
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U  -  «v2)  ■  tE  -  IE  '  E=)  kjfVBj  -  VBE)J  »  E=1  kdv2J  (1-1) 

Is  the  Taylor  expansion  for  the  base-emitter  nonlinear  current  about 
the  operating  point  (I^,Vg^).  The  currents  i^,  1^,  and  ig  are  the 
total  instantaneous,  dc  and  incremental  emitter  currents  and  the 
voltages  vBE>  VBE,  and  v2  are  total  instantaneous,  dc,  and  incremental 
base-emitter  voltages  respectively.  The  nonlinear  current  source  re¬ 
presenting  the  collector  capacitance  nonlinearity  is  given  by 

n 


Yc(vv2)  *  *=1  cj  at  {V3-V2)J  • 


(1-2) 


where  C-'s  are  the  coefficients  of  the  series.  The  collector  dependent 

J 

current  source  is  expressed  as 

g(v2.v3-»2)  •  E__o  9jltv2J(v3-»2)k  .  (1-3) 

k*0 

where  the  g^'s  are  the  coefficients  of  the  series,  and  gCQ=0.  The 
nonlinear  current  source  representing  the  emitter  capacitance  non¬ 
linearity  is 


dv. 


W  =  tF!E  z.s]  W  ^  dt  +  Cje  dt  ’ 


where  is  the  emitter  capacitance  time  constant,  C.  the  base-emitter 

r 

junction  space  charge  capacitance  and  (kT/q)  the  equivalent  thermal 


potential.  The  BJT  device  nonlinearities  therefore  are  represented  by 
series  expansions.  An  example  of  the  series  coefficients  appear  in 


TABLE  1-1 


NONLINEAR  COEFFICIENTS  FOR  2N5L09  BJTa 


Honlineer  Incremental 


Current  Generator 

Coefficient 

Value 

K^Cnfao) 

1.9724 

*  (*2) 

KjCnho/V) 

38.419 

KjCaho/V2) 

498.87 

Cx  (F) 

0.26549  x  10 

*c<VV 

C2  (T/V) 

-0.75664  x  10 

C3  a/v2) 

0.64819  x  10 

«01  <llh0) 

0.22351  x  10 

*10  <aho) 

1.9446 

Cq2 

0. 11426  x  10' 

0 

*20  («ho/V) 

37.754 

,(W-V 

gu  <»ho/V) 

-0.28734  x  10' 

• 

*Q3  (-ho/V2) 

0.30475  x  10 

gjjj  (nho/V2) 

487.88 

gjj  (nho/V2) 

-0.14367  x  10' 

g21  (nho/V2) 

-0.2718  x  10 

Values  reported  by  Fang,  Ref.  18,  p.46. 


Table  1-1.  The  computation  of  the  nonlinear  coefficients  indicated  was 
made  by  the  Nonlinear  Circuit  Analysis  Program  (NCAP)l-^  which  has  a 
built-in  device  model  for  the  BJT.  The  computer  code  requires  only  that 
the  BJT  model  parameters  be  entered.  An  example  of  NCAP  input  para¬ 
meters  for  a  BJT  is  shown  in  Table  1-2.  The  parameters  are  determined 
experimentally^*^  utilizing  several  parameter  extraction  methods. 

The  purpose  of  an  NCAP  analysis  is  to  predict  quantitative  EMI 
effects.  By  comparing  the  NCAP  predicted  results  with  laboratory 
measured  interference  effects,  the  modeling  accuracy  can  be  evaluated. 

An  experimental  system  for  measuring  a  second-order  intermodulation 
effect  is  shown  in  Figure  1-3^^  This  experimental  system  which  simul¬ 
ates  directly  a  practical  interference  situation  has  the  advantage  of 
being  controlled  and  therefore  amenable  to  analysis.  The  experimental 
data  may  be  plotted  as  shown  in  Figure  1-4.  The  data  points  represent 
the  strength  of  the  interfering  signal  at  the  new  frequency  f^  -  f^  as 
a  function  of  the  input  RF  power  at  RF  frequencies  50  kHz,  1  MHz,  and 
10  MHz.  The  AF  voltage  which  is  the  interfering  signal  is  proportional 
to  the  available  input  RF  power,  as  expected,  at  low  RF  power  levels. 

Measurements  can  be  carried  out  at  a  fixed  input  RF  power  while 
the  RF  frequency  is  varied.  The  result  of  such  an  experiment  is  illus¬ 
trated  in  Figure  1-5.  The  most  susceptible  RF  frequency  interval  where 
the  largest  magnitude  interference  is  obtained  is  apparent. 

It  is  desirable  that  the  interference  phenomena  be  expressed  as  a 


TABLE  1-2 


NCAP  INPUT  PARAMETERS  FOR  2N5109  BJTa 


htutttr  Km 


Description 


Talus 


tcb(” 

VCBO  W 
Mu 

Ic  (■*> 

W  (BA) 


1W 

k  (pF-V1). 

b( 

CJ.  (PF, 

C'j  <pF/mA) 

rb  (0) 

rc  (kfl) 

Cj  (pF) 

C3  (pF) 


Avalanche  exponent 
Collector  base  bias  voltage 
Avalanche  voltage 
Collector  capacitance  exponent 
Collector  bias  current 
Collector  currant  at  maxi nun 
Dc  current  gain 
hj.g  nonlinearity  coefficient 
Maximum  d.c.  currant  gain 
Collector  capacitor  scale  factor 
Diode  nonlinearity  factor 
Base-emitter  junction  space 
charge  capacitance 
Derivative  of  base-emitter 
diffusion  capacitance 
Base  resistance 
Collector  resistance 
Base-emitter  capacitance 
Base-collector  and  overlap 
capacitance 


6.0 

5.0 

40.0 

0.285 

50.0 

18.0 

0.363 

84.6 

0.42 

1.0 

10.0 


3.6 


15 

32.9 

0.1 

0.56 


aParameter  values  reported  by  Fang,  Ref.  ig,  p.  45. 
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parameter  of  the  electronic  circuit  or  system  in  question  in  a  manner 
similar  to  the  linear  voltage  gain  parameter.  The  "Nonlinear  Transfer 

r?n 

Function"  derived  from  the  Volterra  Series  kernel1-  J  has  been  used  by 
Graham  and  Ehrman^17^  as  such  a  parameter.  Consider  a  weakly  nonlinar 
network  with  input  x(t)  and  output  y(t).  Let  X(f)  and  Y(f)  be  the 

Fourier  transforms  of  x(t)  and  y(t).  Then  Y(f^  is  related  to  X(f)  by  the 

•  [17] 
expression1-  J 

+a> 

Y(f)  -  Z  JJ...J  Hn(fr...,fn)  X(f1)...X(fn) 


•«  (f-frf2"--*"fn)  dfl  df2---dfn 


(1-5) 


In  the  above  equation  Y(f)  is  the  sum  of  all  orders  of  nonlinear  re¬ 
sponses;  X(f)  is  the  input  signal  spectrum,  6(f-f1-f2~. . .-fn)  is  the 
delta  function  defined  by  the  integral 

a> 

«(f-frf2-...-fn)  =  f  exp  {-j2ir(f~frf2-...-fn)t}  dt  .  (1-6) 

*  -• 

and  Mn(f|,f2»-.*»fn)  is  the  nth  order  nonlinear  transfer  function.  The 
first  order  (linear)  transfer  function  can  be  expressed  as^7^ 


H*| (f i )  =  [Y (f i )]"  JA(f.),  i  -  1,2.. ...n 


(1-7) 


where  Y(f^)  is  the  network  admittance  matrix  evaluated  at  f . ,  and 
JA(f i )  the  first  order  current  source  vector  that  excites  the  network. 
The  second  order  nonlinear  transfer  function  H2(f^,f2)  can  be  expressed 

..nn 


WV  ■ <1/2)  [y(Vfj)]  Jd(W  ; 

for  i  =  1,2,. ...n;  j  =  1,2,. ..,n 


(1-8) 


.*•  :•  %’  \m  •  •  *v\  v,.-. v, . 


where  Y (f . )  Is  the  network  admittance  matrix  evaluated  at  the 
frequency  (f.+f.),  and  Jd(f<.f.j)  is  the  second  order  nonlinear  current 

I  J  si  l  J 

source  vector  that  excites  the  network.  The  third  order  nonlinear 
transfer  function  H3(f^  ,‘f »f k)  can  be  expressed  as^18^ 

H3(fi'fj’fk}  =  $  MVWJ"1  Jd(fi  *fj*fk); 

1  *  1,2,. ...n;  j  -  1,2, ...,n;  k  =  l,2,...,n  (1-9) 


where  Y(f.+f .+f. )  is  the  admittance  matrix  evaluated  at  the  frequency 

I  J  lx 

(f  1+f j+f k) *  and  Jjlf-j ,fj»fkJ  is  ^e  avara9a*  nonlinear  third  order  cur¬ 
rent  source  vector  that  excites  the  network.  Higher  order  nonlinear 
transfer  functions  can  be  expressed  in  a  similar  manner. 


The  experimental  system  of  Figure  1-3  is  used  to  measure  the 
second  order  nonlinear  transfer  function.  Let  the  two  input  signal 
frequencies  f-j  and  f2  be  denoted  as  f2  =  fR-p  and  f^  =  (fRp  -  fAp)  where 
fRp  Is  the  RF  carrier  frequency  and  f^p  is  the  AM  modulation  frequency. 


The  AF  voltmeter  is  tuned  to  the  frequency  fAp  =  f2  -  =  f2  +  (-f1 ) . 

The  component  of  the  output  voltage  at  f2  -  f-j  can  be  expressed  as^8^ 


v(t)  «  mA^  H2(f2,  -  fj)  cos[2Tr(f2  -  f^t] 


(1-10) 


where  A  is  the  amplitude  of  the  RF  carrier  at  f2  =  fRp,  and  m  is  the 
modulation  index.  The  rms  voltage  read  on  the  AF  voltmeter  is 


V  =  0.707mA*  iH^.-f^l 


(1-11) 


* 


For  the  third  order  current  vector  averaging  process,  see  Chapter  5 
of  Reference  15. 


Expressed  in  dB  with  respect  to  a  1  mV  reference  level.  Equation  1-11 
becomes 

201og|v/lmV|  =  201og[0.707  mA^H^fg.-fj)  |]  -  201og  (10“3) 

«  57  +  201  ogm  +  401ogA  +  201og |H2(f2,-f-, }  j  (1-12) 

When  the  RF  generator  impedance  equals  50  n,  the  RF  carrier  amplitude 
A  is  related  to  the  available  generator  power  by  P  =  A2/400.^13^ 

With  m  =  0.5,  Equation  (1-12)  can  be  written  as 

201og|V/lmV|  =  2Pgen(dbm)  +  43  +  201og|H2(f2,-f1) |  ,  (1-13) 

or 

201og|H2(f2,-f1)|  =  201og|V/lmV|  -43  -  2Pgen(dBm)  (1-14) 

Equation  1-13  indicates  that  the  AF  voltmeter  reading  V  is  proportional 
to  Pgen.  At  low  values  of  Pgen  the  data  plotted  in  Fig.  1-4  obey 
Equation  (1-13).  Using  Equation  (1-14)  and  the  data  in  Figure  1-4  at 
low  Pgen  values,  values  of  the  second  order  transfer  function  H2(f2,-f1) 
can  be  determined.  Values  of  H2(f2>-f^}  vs  frequency  f2  are  plotted  as 
illustrated  in  Figure  1-6.  It  is  to  be  noted  that  in  both  Figures  1-5 
and  1-6  the  frequency  (-fj)  is  varied  in  step  with  f2  to  keep  (f2  -  f ^ ) 
a  constant  at  a  value  fAF  ’  f2  -  fr 

An  experimental  system  for  measuring  a  third  order  nonlinear 
transfer  function  is  shown  In  Figure  1-7. ^20-^  For  an  electronic  circuit 
with  third  degree  nonlinearities,  two  tone  excitation  can  give  rise  to 
the  following  third  order  frequency  mixes:  (2^-^),  (2f2-f1),  (3^), 
(2f,+fp),  (2f?+f.),  and  (3f?).  Since  any  one  of  the  third  order 


frequency  mixes  is  representative  of  the  circuit  third  degree  nonlinear¬ 


ities,  only  the  output  at  one  such  frequency  2f2  -  f-|  will  be  considered 
as  an  example.  The  spectrum  analyzer  rms  voltage  V( 2f^  -  f^)  at  fre¬ 
quency  2fg  -  fj  is  related  to  the  third  order  nonlinear  transfer  func¬ 
tion  H3(-f^  .^.fg)  by  the  expression^1 8^ 

V(2f2  -  f,)  -  (3/4)|H3(-frf2,f2)|A1A22  ,  (1-15) 

where  A^  and  A2  are  the  voltage  amplitude  of  the  two  input  signals 
at  frequencies  f^  and  f,,. 

1.2  System  Nonlinearity  of  the  Integrated  BJT  Operational  Amplifier 

The  application  of  integrated  circuits,  especially  the  integrated 
BOT  Op  Amp,  permeates  modern  electronic  equipment  design.  The  EMI 
analysis  of  the  IC  Op  Amp  has  become  an  important  consideration.  Such 
an  analysis  poses  a  number  of  challenges.  The  first  is  the  complexity 
of  the  circuitry  where  20  to  30  BJT's  are  interconnected.  The  second 
Is  the  need  to  know  the  nonlinear  parameters  of  each  and  every  indivi¬ 
dual  BOT.  Finally,  the  effects  of  the  parasitic  elements  associated 
with  the  BOT's  have  to  be  accounted  for.  Recently  the  nonlinear 
analysis  of  the  uA741  Op  Amp  was  successfully  achieved  by  Fang.^18^ 

The  schematic  circuit  diagram  of  this  widely  used  IC  is  shown  in 
Figure  l-8.^23^  The  Nonlinear  Circuit  Analysis  Program  (NCAP)  was  the 
main  analysis  tool  used  by  Fang.  The  NCAP  input  coding  diagram  for 
the  741  buffer  amplifier  is  shown  in  Figure  1-9. ^18^  The  NCAP  pre- 
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Figure  1-8.  Schematic  Diagram  of  the  7^1  Op  Amp 


diction  for  the  first  order  transfer  function  vout/VIN  vs  frequency  is 
compared  with  the  experimental  data  in  Figure  1-10.  The  NCAP  calcula¬ 
tion  of  the  second  order  nonlinear  transfer  function  ^(fg.-f^)  vs  RF 
frequency  is  compared  with  the  experimental  data  in  Figure  1-11.  It 
can  be  seen  that  the  nonlinear  analysis  utilizing  NCAP  predicts  accurate¬ 
ly  the  location  of  the  most  susceptible  RFI  RF  frequency  range  and  the 
general  shape  of  the  RFI  frequency  response.  The  accuracy  achieved 
establishes  the  nonlinear  analysis  techniques  as  a  useful  tool  for  EMI 
prediction. 

The  NCAP  input  coding  diagram  shown  in  Figure  1-9  with  104  nodes 
creates  a  very  large  admittance  matrix.  The  NCAP  computer  code  employs  a 
sparse  matrix  technique  to  reduce  the  in-core  storage.  The  data  are  pass¬ 
ed  back  and  forth  between  the  Central  Processing  Unit  (CPU)  and  the  pe¬ 
ripheral  accessory  memories.  For  a  five-frequency  sweep,  determining  the 
741  buffer  amplifier  second  order  transfer  function  requires  approximately 
half  an  hour  CPU  time  and  another  half  hour  peripheral  time  on  a  Honeywell 
6180  computer.  Up  to  500  nodes  can  be  handled  by  the  present  NCAP 
algorithm.  The  computational  procedure  used  by  NCAP  begins  with  the  con¬ 
struction  of  nodal  equations  using  Kirchhoff's  Voltage  Law.  The  first 
order  transfer  function  in  Equation  (1-7)  is  calculated  next  at  each 
frequency  f^  required.  The  second  order  nonlinear  current  source  vector 
Jd(^i»f|c)  is  then  determined  at  the  frequency  (f^  +  f^)  using  the  first 
order  nodal  voltages  at  frequencies  f^  and  f^.  The  appropriate  second 


Second  Order  Transfer  Function  H^rfgXdB)  3  VOUT/VIN  (dB) 


order  transfer  function  is  then  calculated  from  Equation  (1-8).  The 
procedure  is  repeated  for  each  combination  of  frequencies  (f.,f .)  re- 

*  w 

quired.  Similarly,  the  third  and  higher  order  transfer  functions  are 
calculated.  The  NCAP  program  can  calculate  up  to  the  sixth  order  trans¬ 
fer  function.  It  is  to  be  noted  that  a  calculation  of  a  third  order 
transfer  function  for  a  two  tone  excitation  requires  that  first  order 
transfer  functions  at  f^  and  fg  and  second  order  transfer  functions  at 
fg  -  fj,  2f-|,  f-j  +  f^  and  2fg  be  calculated  first. 

The  NCAP  model  parameters  for  the  BJTs  in  the  pA741  were  deter¬ 
mined  by  Fang^®]  in  two  ways.  Fang  used  the  modified  Ebers-Moll  BJT 
model  parameters  reported  by  Wooley  et  al.^^  for  the  Intersil  741 
type  chip.  Fang  also  measured  parameters  in  situ  on  a  Fairchild 
pA741  devices  using  probe  methods. Input  parameters  for  both  the 
SPICE2  (Simulation  Program  with  Integrated  Circuit  Emphasis)^53-! 

Computer  Program  and  NCAP  were  determined.  There  are  six  BJT  types  in 
the  741  IC:  small  npn,  large  npn,  lateral  pnp,  large  substrate  pnp, 
dual -col lector  lateral  pnp  and  dual -emitter  substrate  pnp.  The  "small" 
and  "large"  refer  to  the  device  geometry,  ranging  from  60  pm  x  100  pm 
to  140  pm  x  200  pm.  For  the  first  time  the  NCAP  BJT  model  parameters 
for  Integrated  BJTs  became  available.  The  application  of  nonlinear 
analysis  computer  program  NCAP  for  EMI  analyses  of  bipolar  ICs  was  made 
possible.  This  was  a  major  achievement.  Fang's  NCAP  input  coding  data 
are  reproduced  in  Table  1-3.  Fang  reported  that  there  is  considerable 
uncertainty  in  the  measurement  of  substrate  capacitances  in  both  his  own 


and  Wooley's  data.  The  individual  substrate  capacitance  values  may  not 
correspond  directly  to  the  device  surface  geometry.  The  important 
effects  of  the  substrate  capacitance  values  on  the  nonlinear  analysis 
of  IC  Op  Amps  are  considered  in  detail  in  this  research. 

1.3  Macromodeling  of  BJT  ICs 

The  ICs  fabricated  today  are  complex,  and  the  ICs  fabricated  to¬ 
morrow  will  be  more  complex.  The  complexity  of  ICs  makes  computer-aided 
analysis  of  EMI  in  IC  circuits  time  consuming  and  expensive.  Further-, 
more,  the  EMI  analyst's  comprehension  of  the  overall  system  behavior 
tends  to  be  obscured  when  he  attempts  to  keep  track  of  what  is  happening 
at  the  device  level.  There  is  a  need  of  a  simplified  and  yet  re¬ 
presentative  equivalent  circuit  for  the  1C.  This  need  has  lead  to 

the  development  of  macromodels. 

There  are  a  number  of  approaches  to  macromodeling.  The  IC  may  be 
represented  by  an  approximation  of  one  or  more  of  its  transfer  functions 
Another  method  is  to  represent  the  IC  by  its  terminal  behavior.  The 
goal  of  macromodeling  is  to  simplify  the  internal  IC  circuitry.  The 
result  can  be  reduction  of  computation  time  required  for  an  analysis  by 
four  to  ten  fold.  When  IC  macromodels  are  used,  a  system  employing  many 
IC's  becomes  amenable  to  analysis  by  the  circuit  designer  and  the  EMI 
analyst.  Macromodeling  efforts  all  share  some  common  characteristics. 
The  range  of  application  is  necessarily  limited.  Nearly  all  macromodels 


are  developed  in  conjunction  with  computer  aided  circuit  design  tech¬ 
niques  and  are  intended  for  computer  analysis.  One  clear  distinction 
need  be  made  in  our  definition  of  a  macromodel.  That  is,  a  macromodel 
Involves  an  approximate  representation  of  one  or  more  actual  transistors 
in  an  IC.  Therefore  an  idealized  transfer  function  representation, 
such  as  a  voltage  gain  parameter  to  represent  an  amplifier  for  example, 
is  not  a  macromodel  in  the  context  of  the  present  research. 

Narud  reported  a  macromodel  for  digital  ICs  in  1970. This 
early  model  maintained  the  individual  device  characteristics  at  the 
inputs,  up  to  the  gate  or  flip-flop  level.  The  balance  of  the  IC  is 
represented  by  an  idealized  black-box.  This  approach  has  proven  to  be 
widely  useful.  Digital  IC  macromodels  have  been  developed  continuously 
and  now  are  widely  used  in  the  analysis  and  design  of  large  scale 
integration  (LSI)  IC's.^»28]  However,  the  various  digital  macromodels 
used  by  different  manufacturers  are  usually  proprietary  information 
and  are  seldom  available  to  the  EMC  community. 

Boyle  et  al.  reported  a  macromodel  for  the  IC  Op  Amp  which  can  be 
used  for  nonlinear  dc  and  transient  analysis. This  macromodel  pro¬ 
vides  a  pin-for-pin  representation  as  shown  in  Figure  1-12.  The  values 
of  the  elements  in  the  macromodel  are  determined  from  the  Op  Amp 
specifications  and  input-output  responses.  The  three-stage  model  in¬ 
clude  elements  related  to  dc  circuit  quantities:  I£E  for  the  dc  cur¬ 
rent  source  in  the  input  stage,  R£  related  to  the  Early  Voltage  of 
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Figure  1-12.  A  Macroraodel  for  the  7^1  Type  BJT  IC 

Op  Amp  for  DC,  Small-Signal  and  Trans¬ 
ient  Analyses  Developed  by  Boyle  et  al 


the  current  source  transistor,  Rp  representing  the  dc  power  dis¬ 
sipation,  output  stage  current  limiting  parameters,  ,  R^,  D-j ,  D2  and 
RqP  and  the  output  voltage  limiting  elements  V^,  Dg,  V^,  and  D^.  Small - 
signal  operation  is  modeled  by  the  elements  G_  for  the  differential 

a 

gain,  Gcm  for  the  Common  Mode  Rejection  Ratio  (CMRR),  G^  for  the  inter¬ 
stage  gain,  Rgg  for  the  ac  output  resistance,  for  the  differential 
amplifier  excess  phase  effect,  and  Cg  for  the  3  dB  cutoff  frequency. 

The  positive  slew  rate  is  related  to  Cj»  and  the  negative  slew  rate  to 
and  C2- 

It  is  significant  that  Boyle  et  al.  found  that  four  junctions  at 
the  input  stage  are  necessary  and  sufficient  to  model  the  741  IC  per¬ 
formance.  This  is  in  agreement  with 'the  work  of  Narud.  The  Boyle 
macromodel  employs  two  ideal  transistors  which  contain  four  p-n  junctions 
An  earlier  model  by  Treleavan  and  Trofimenkoff ,^30^  shown  in  Figure  1-13, 
employs  no  representation  of  any  physical  junctions.  Other  macro- 
models^-33^  also  eschew  the  direct  inclusion  of  physical  devices  as 
shown  in  Figure  1-14.  The  necessity  for  physical  modeling  in  the  non¬ 
linear  application  of  the  macromodel  will  be  developed  in  this  disserta¬ 
tion. 

The  macromodel  developed  by  Boyle  et  al.  can  be  used  for  dc,  small- 
signal,  and  transient  calculations.  However,  it  was  not  designed  for 
nonlinear  analysis  related  to  EMI  predictions.  A  macromodel  has  also 
been  developed  for  the  IC  Comparator  J 34 Usually,  about  an  order  of 


magnitude  reduction  In  matrix  size  and  computation  time  is  realized 
when  these  macromodels  are  used  in  place  of  device  level  or  full 
models. 

For  completeness,  two  other  nonlinear  macromodels  for  the  Op  Amp 
should  be  mentioned  even  though  no  attempt  has  been  made  to  use  these 
nonlinear  macromodels  for  EMI  analysis.  The  Weil  nonlinear  macro- 
model  represents  the  nonlinearities  by  two  diodes  at  the  input  of 
the  IC.  The  diode  current-voltage  characteristic  is  exponential  and 
includes  an  ideality  factor.  The  input  resistance  is  nonlinear  as  a 
result  of  the  diode  exponential  current-voltage  characteristic.  The 
nonlinear  macromodel  developed  by  Majewski  and  Sanchez-Sinencio^-3®’37^ 
falls  in  the  category  of  transfer  function  modeling  in  the  manner  of 
piecewise  linear  analysis, f39,40^  as  opposed  to  physical  device  modeling. 
In  their  model,  the  saturation  or  clipping  of  the  input  current  versus 
input  voltage  represents  the  overall  IC  nonlinearity.  Output  waveform 
distortion  under  large  signal  operation  can  be  predicted  with  the  aid 
of  Fourier  analysis.  Output  voltage  limiting  is  represented  by  a  delta 
function  type  of  nonlinear  resistance.  These  nonlinear  macromodels  were 
not  designed  for  EMI  applications  nor  can  the  BJT  nonlinear  current 
generators  shown  in  Figure  1-2  be  included  in  these  models. 

Macromodels  for  ICs  have  been  reported  upon  by  several  investi¬ 
gators.^"3®^  The  macromodel  for  the  Op  Amp  developed  by  Boyle  et  al. 
was  selected  for  EMI  analysis  because  it  includes  two  transistors  which 


model  In  a  physical  sense  the  actual  Op  Amp  input  circuitry.  Because 
the  Boyle  Op  Amp  macromodel  predicts  EMI  in  real  Op  Amps  so  well,  no 
effort  to  develop  or  to  modify  the  other  IC  macromodels  for  EMI  analyses 
was  made.  The  remainder  of  this  dissertation  is  devoted  to  describing 
how  well  the  Boyle  Op  Amp  macrociodel  predicts  the  demodulation  of  AM 
modulated  RF  carriers  in  op  amp  circuits  and  to  explaining  why  it  works 
so  well. 

The  results  of  the  research  reported  in  this  dissertation  has  been 
summarized  in  a  paper  which  was  presented  at  the  International  Conference 
on  Electromagnetic  Compatibility  at  Southhampton,  England,  the  16th  to  18th, 
September,  1980.  A  copy  of  this  paper  is  attached  as  Appendix  A. 


CHAPTER  TWO 


THE  NONLINEAR  MACROMODEL  FOR  EMI  ANALYSIS  OF  THE 
BIPOLAR  IC  OPERATIONAL  AMPLIFIER 

As  discussed  in  Chapter  1  the  Op  Amp  macromodel  developed  by 
Boyle  et  al.  was  selected  for  EMI  analysis  because  it  includes  two 
transistors  which  model  in  a  physical  sense  the  Op  Amp  input  circuitry. 
In  Section  2.1  the  transformation  of  the  original  Boyle  macromodel, 
which  was  developed  primarily  for  dc  and  transient  nonlinear  analyses. 
Into  a  nonlinear  small -signal  macromodel  suitable  for  use  with  NCAP 
will  be  described.  In  Section  2.2,  intermodulation  predictions  made 
using  the  nonlinear  small -signal  Op  Amp  macromodel  will  be  compared 
to  predictions  made  using  a  full  Op  Amp  model  and  to  experimental 
results.  The  specific  case  investigated  involves  the  demodulation  of 
amplitude  modulated  RF  signals  in  Op  Amps  to  produce  undesired  low 
frequency  responses  related  to  AM  modulation  on  the  RF  carrier.  In 
Section  2.3  the  important  effects  caused  by  parasitic  capacitances 
associated  with  the  p-n  junctions  used  to  provide  isolation  in  in¬ 
tegrated  circuit  Op  Amp  is  discussed  in  some  detail. 

2.1  The  Nonlinear  Macromodel  Based  on  a  Two-Transistor  Physical 

Representation 

An  analog  IC  can  be  blocked  into  three  main  sections:  ^  the 
input  stage,  the  gain  stage  and  the  output  stage.  For  the  purpose  of 


nonlinear  small -signal  analysis,  the  output  stage  can  be  viewed  as  a 
linear  stage  that  provides  power  amplification  and  impedance  trans¬ 
formation.  The  input  stage  is  a  differential  amplifier.  The  gain 
stage  provides  nearly  all  of  the  open-loop  gain  and  phase  compensation. 

A  schematic  representation  of  the  IC  Op  Amp  can  be  drawn  as  in  Figure 
2-1. As  noted  previously  this  two- transistor  configuration  fur¬ 
nishes  the  means  of  including  physical  device  nonlinearities  because 
the  BJT  nonlinear  current  generators  shown  in  Figure  1-2  can  be  included. 
The  Op  Amp  model  shown  in  Figure  2-1  is  a  simple  model.  It  can  be  re¬ 
fined  as  necessary.  It  is  interesting  to  note  that  the  use  of  a  two-  . 
transistor  input  configuration  followed  with  a  linear  amplifier  for  an 
analog  IC  was  used  in  an  early  computer  simulation. 

In  the  nonlinear  small-signal  macromodel  shown  in  Figure  2-2,  the 
transistors  Q1  and  Q2  are  the  two  BJTs  at  the  inputs  of  the  IC  Op  Amp 
being  modeled.  These  two  BJTs  dan  be  characterized  by  the  nonlinear 
BJT  model  of  Figure  1-2  and  by  the  model  parameters  of  Table  1-3.  The 
nonlinear  behavior  of  the  Op  Amp  is  caused  by  these  two  transistors. 

The  nonlinear  macromodel  of  Figure  2-2  can  be  derived  directly  from  the 
Boyle  macromodel ^3  of  Figure  1-12.  When  the  nonlinear  small-signal 
behavior  of  the  Op  Amp  is  to  be  modeled,  several  dc  and  limiting  circuit 
elements  in  the  Boyle  macrocode!,  namely  Rp,  Dp  D2,  Rc,  Gcvg,  Vc,  D3, 
and  Vg  can  be  omitted.  In  the  small-signal  macromodel  the  capacitor 
C-j  is  related  to  the  nondominate  pole  of  the  differential  gain  response 


Figure  2-1.  Schematic  Representation  of  the  741  Type  Op  Amp. 

The  signal  generator  F  is  a  function  of  the  input 
voltage  V^. 


Figure  2-2.  Nonlinear  Macromodel  for  the  BJT'  IC  Op  Amp 


in  Small-Signal  Operation. 
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and  to  the  excess  phase  at  the  cutoff  frequency  fQdg.  The  capacitor  C2 
is  the  phase  compensation  capacitance  in  the  g?.in  stage  of  the  actual 
IC.  The  resistors  =  Rc2  are  determined  by  the  requirement  that  the 
0  dB  frequency  f0dB  =  l/2itRcC2.  The  resistors  Rel  "  Re2  model  the  emit¬ 
ter  degeneration  for  slew  rate  enhancement. [47,48]  y^g  res-jstor  R£  is 
related  to  the  differential  input  resistance.  The  capacitor  models 
charge  storage  effects  at  the  input.  The  transconductances  Ga  and  Gfa 
are  related  to  the  differential  voltage  gain.  The  value  for  the  volt¬ 
age  gain  of  the  differential  input  stage  is  not  critical  in  either 
linear  or  nonlinear  analysis  and  is  assigned  a  value  of  unity.  The  in¬ 
put  stage  differential  mode  voltage  gain  equals  the  product  G  x  R0. 

The  product  G^  x  Rq2  equals  the  single-ended  open-loop  differential 
'  voltage  gain  of  the  IC.  The  resistor  RQj  is  the  dc  component  of  the 
internal  output  resistance.  The  resistor  RQ2  is  a  small -signal  inter¬ 
nal  output  resistance.  The  common-mode  transconductance  G  models 

cm 

the  effect  of  the  Common-Mode  Rejection  Ratio  (CMRR). 

* 

The  macromodel  circuit  parameters  can  be  extracted  from  the  IC 
manufacturer's  specification  sheets  using  procedures  given  by  Boyle  et 
al.  The  macromodel  parameters  for  the  LM741  Op  Amp  are  given  in  Table 
2-1. 


The  NCAP  BJT  model  parameters  for  the  input  transistors  Q1  and  Q2 
In  the  macromodel  are  taken  directly  from  the  published  data  by  Fang 
and  Whalen  and  are  reproduced  in  Table  2-2.  The  BJT  parameters  used 
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TABLE  2-1 

MACROMODEL  PARAMETERS  FOR  LM741a 


300®  K 

Ga 

188.6  umho 

8  *  10'16A 

GCM 

6.23  nnito 

100  kC 

«oi 

32.13  a 

30  pF 

*02 

42.87  Q 

2.41  pF 

Sb 

247.49  mho 

20.26  uA 

1SD1 

8  x  10‘16A 

9.872  m 

*C 

0.02129  x  10' 

5305  Q 

Gc 

49964  mho 

2712  Q 

vc 

1.803  V 

5.460  pF 

- 

VE  * 

2.303  V 

Erora  Boyle  et  al.,  Ref.  29 


TABLE  2-2 

NCAP  INPUT  PARAMETER  VALUES  FOR  pA?4l  OP  AMP  MACROMQDEL3, 


Parameter 

Q1 

Q2 

n 

4.34 

4.34 

VCB  ^ 

14.57 

14.57 

vCB0  ^ 

20 

20 

U 

0.165 

% 

0.165 

k  (A) 

10,E-6 

10.E-6 

*C  max  W) 

0.5. E-3 

0.5  E-3 

a 

0.90 

0.62 

hFE  max 

400 

400 

k  (F-  V  \) 

1.23  E-12 

1.23  E-12 

Ref. 

1.091 

1.091 

Cje  (F) 

T.23  E-12 

1.23  E-12 

(F/A) 

9.09  E-9 

9.09  E-9 

rb  (ohms) 

830 

830 

rc  (ohms) 

5.33  E6 

5.33  E6 

(F) 

0.1  E-12 

0.1  E-12 

C3  (F) 

0.1  E-12  • 

0.1  E-12 

*This  set  of  parameter  values  is  'based  upon  the  full  Op  Amp 
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BJT  parameter  values  reported  by  Fbng,  with  the  exceptions 
of  I  and  a  as  noted  on  page  40. 


are  those  for  the  Op  Amp  input  transistors  which  are  small  npn  transis¬ 
tors.  As  will  be  discussed  later,  the  nonlinear  BJT  parameters  of 
Table  2-2  can  be  considered  applicable  to  BJT  IC  Op  Amps  of  present 
day  technology. 

Utilizing  the  nonlinear  small-signal  macromodel  described  in  this 
section,  the  NCAP  predictions  for  the  second  order  transfer  function 
of  the  741  unity  gain  buffer  amplifier  are  reported  upon  in  the  next 
section. 

2.2  NCAP  Intermodulation  Predictions  with  the  Nonlinear  Macromodel  and 
Experimental  Verification. 

* 

The -validity  of  the  nonlinear  macromodel  proposed  in  Figure  2-2 
has  to  be  established  by  determining  its  effectiveness  in  EMI  analysis. 
The  nonlinear  circuit  analysis  program  NCAP  is  used  for  the  validation. 
Two  criteria  for  the  nonlinear  macromodel  are  proposed: 

1.  Do  NCAP  predictions  based  upon  the  macromodel  agree  with 
the  same  predictions  based  upon  the  full  model  (device 
level  model)  of  the  Op  Amp  IC? 

2.  Do  NCAP  predictions  of  EMI  effects  based  upon  the  macro¬ 
model  agree  with  the  experimental  results? 

The  741  Op  Amp  circuit  selected  for  the  validation  is  the  unity  gain 

buffer  amplifier.  The  circuit  configuration  is  shown  in  Figure  2-3.  T 

nonlinear  macromodel  equivalent  circuit  for  this  amplifier  is  shown  in 

Figure  2-4.  The  experimental  system  for  the  measurement  of  the  second 

order  transfer  function  H?(f?,-f.)  is  shown  in  Figure  2-5. 


The  coding  diagram  for  an  NCAP  analysis  with  the  nonlinear  macro- 
model  is  shown  in  Figure  2-6.  The  circuit  of  Figure  2-6  represents  the 


small-signal  equivalent  circuit  of  the  experimental  circuit  of  Figure 
2-5.  If  nonlinear  effects  do  not  occur  in  the  IC  Op  Amp,  the  Wave 
Analyzer  shown  in  Figure  2-5  should  read  zero.  Then  the  NCAP  EMI 
predictions  for  Figure  2-6  should  be  zero  also.  It  is  noted  that  any 
EMI  predicted  for  the  macromodel  circuit  of  Figure  2-6  results  solely 
from  the  nonlinearities  of  the  two  BJT's.  The  NCAP  input  coding  list 
for  a  second  order  intermodulation  analysis^^  of  the  circuit  of  Figure 
2-6  is  given  in  Table  2-3.  The  circuit  element  parameter  values  used 
In  Table  2-3  are  taken  from  Figure  2-6.  The  nonlinear  BJT  parameter 
values  given  in  Table  2-2  are  taken  from  Table  1-3  with  two  exceptions. 
Fang^^  reports  a  value  1^  =  7.85  jiA.  Boyle  et  al.  uses  a  value 
Ic  *  10  uA.  The  Boyle  value  is  used  in  our  analysis.  The  value  of  the 
hpE  Nonlinearity  Coefficient  "a"  is  adjusted  using  the  following 
empirical  relationship: 


a 


hFEmax. ~ 

I 

hpE[log(j 


Cmax 


)] 


(2-1) 


The  values  obtained  for  "a"  of  the  macromodel  transistors  and 

are 


ai 


400  -  111.67 


0.9 
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111.67  [log ( 10,0 

0.5  x  10" 


Figure  2-6.  NCAP  Coding  Circuit  Diagram  for  the 
Voltage  Follower  Circuit  with  Op  Amp 
Macromodel . 


The  0.01  0  resistors  represent 
jumper  wires,  as  0  0  resistors  are 


not  allowed  in  NCAP. 


Table  2-3 

NCAP  INPUT  CODE  FOR  SECOND  ORDER  INTERMODULATION  ANALYSIS 
OF  THE  EXPERIMENT  OF  FIGURE  2-5 


**  MACROMODEL  OF  741  OP  AMP 
**  ANALYSED  BY  NCAP  FOR  H2 
+*  SUBMITTED  AT  SUNY  BUFFALO 
**  GORDH  OUTPUT  NODE  IS  15 

*  START  CIRCUIT 

*  GENERATOR 
NODE  1  0 

FR  1  1E6  10E6  5  LIN 
AMP  1.0  0.0 

FR  2  -0.9996E6  -9.9996E6  5  LIN 
AMP  1.0  0.0 
IMP  50  0 

*  LINEAR  COMPONENTS 
R  2  1  0.01 

R  2  3  560 
R  3  8  0.01 
R  2  0  50 
R  0  6  5305 
R  7  16  2712 
R  0  10  5305 
R  11  16  2712 
R  16  0  9.872E6 
R  12  0  1.0E5 
R  13  jO  42.87 
R  13  14  32.13 
R  14  15  4.7E4 
R  15  0  4.489E4 
R  14  4  620 
C  2  0  1.1E-10 
C  6  10  5.460E-12 
C  16  0  2.41E-12 


C  12  13  3.0E-11 
C  15  0  2.106E-9 
**  SUBSTRATE  CAPACITANCE 
C  6  0  4.0E-1 5+ 

C  10  0  4.0E-15 
C  7  0  4.0E-15 
C  11  0  4.0E-15 

*  TRANSISTOR 
NODE  4 

4.34  14.57  20  0.165  10.05E-6 
0.5E-3  0.90  400  1.23E-12  1.091 
1.23E-12  9.09E-9  830  5.33E6 
0. IE-12  0.1E-12 

*  TRANSISTOR 
NODE  8 

4.34  14.57  20  0.165  10.05E-6 
0.5E-3  0.62  400  1.23E-12  ’ 

1.091  1.23E-12  9.09E-9  830 
5.33E6  0.1E-12  0.1E-12 

*  LINEAR  DEPENDENT  SOURCE 
NODE  6  10  0  12 

VC  1.886E-4  0.0 

*  LINEAR  DEPENDENT  SOURCE 
NODE  16  0  12  0 

VC  6.28E-9  0.0 

*  LINEAR  DEPENDENT  SOURCE 
NODE  12  0  0  13 

VC  247.49  0.0 

*  END  CIRCUIT 

*  END 


t  For  ease  of  adding  and  eliminating  substrate  capacitances  in  NCAP 
runs,  a  0  pF  substrate  capacitance  is  represented  by  a  reduction 
of  three  orders  of  magnitude. 
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=  0.62 
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a2  * 


400  -  143.57 


,-6 


143.57  [log(10>0  x  1  V~)]2 
0.5  x  10‘J 

The  values  for  the  second  order  transfer  function  ^(fp-fg) 

calculated  by  NCAP  as  a  function  of  RF  frequency  are  plotted  in  Figure 

2-7.  In  order  to  generate  the  transfer  function  spectrum  plotted,  the 

frequency  sweeping  feature  of  the  NCAP  program  was  used.  The  experi- 
♦ 

mental  results  obtained  using  the  system  shown  in  Figure  2-5  were 
converted  to  experimental  values  for  the  second  order  nonlinear  transfer 
function  using  Equation  1-14  and  are  plotted  in  Figure  2-7. 

The  results  shown  in  Figure  2-7,  while  not  satisfactory,  are  en¬ 
couraging.  The  macromodel  does  predict  the  frequency  range  at  which 
maximum  second  order  intermodulation  occurs.  Although  the  predicted 
second  order  transfer  function  is  in  error  by  15  dB,  there  is  qualitative 
agreement  between  macromodel  predictions  and  experimental  results. 

Plotted  on  the  same  figure  are  the  NCAP  predictions  for  the  same  buffer 
amplifier  based  on  a  full  model  for  the  Op  Amp.  It  Is  seen  that  the 
full  model  predicted  results  are  higher  than  the  macromodel  by  5  dB  in 
the  critical  frequency  range  of  maximum  second  order  intermodulation 

*  The  experimental  results  were  obtained  using  the  Fairchild  yA741 
IC  which  is  essentially  the  same  as  the  LM741 . 

**  Values  for  a  third  order  nonlinear  transfer  function  were  also  cal¬ 
culated  and  determined  to  be  less  than  -100  dB.  An  experimental 
verification  of  such  a  low  intermodulation  signal  strength  would 
be  very  difficult  if  not  impossible  and  was  not  attempted. 
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effects,  but  are  still  10  dB  lower  than  the  experimental  results.  Similar 
NCAP  predictions  by  Fang  et  al.^®^  using  the  full  model  for  the  Op  Amp 
matches  the  experimental  results  more  closely.  The  Fang  NCAP  simulation 
accounted  for  the  parasitic  effects  of  the  substrate  capacitances  (See 
Figure  1-9).  These  parasitic  effects  are  not  accounted  for  in  the 
NCAP  predictions  of  Figure  2-7.  The  need  to  account  for  the  parasitic 
effects  of  the  substrate  capacitances  in  nonlinear  IC  Op  Amp  modeling 
is  evident.  Substrate  capacitance  effects  are  discussed  in  the  next 
section. 

2.3  Substrate  Capacitances  in  Nonlinear  Macromodeling 

In  the  EMI  modeling  of  electronic  circuits,  Whalen  and  Paludi^49-^ 
showed  that  the  parasitic  elements  strongly  influence  the  accuracy  of 
computer-aided  analysis.  This  is  the  case  even  for  circuits  containing 
only  one  active  device  which  is  the  sole  source  of  circuit  nonlinearity. 

riel 

FangL  J  in  his  computer  simulation  of  Op  Amps  using  a  full  model  in¬ 
cluded  substrate  capacitances  as  parts  of  his  IC  circuit  model.  The 
substrate  capacitances  of  the  741  measured  by  Wooley  et  al.^^  were 
used.  Also,  a  lower  typical  value  of  2  pF  was  used.  Consequently,  Fang 
was  able  to  make  predictions  of  intermodulation  effects  closely  matching 
experimental  results.  Although  Wooley  found  that  the  substrate  capaci¬ 
tances  limits  the  first  order  response  significantly,  none  of  the  macro¬ 
models  referenced  in  Section  1.3  include  these  parasitic  elements. 

Since  those  macromodels  were  not  designed  for  EMI  nonlinear  analysis,  the 


effects  of  neglecting  substrate  capacitances  were  not  appreciated. 

The  npn  transistors  in  a  bipolar  IC  Op  Amp  have  substrate  capac¬ 
itances  occurring  at  the  collector.  These  capacitors  are  shown  in  the 
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741  full  model  NCAP  coding  diagram.  Figure  1-9.  According  to  Wooley,L  J 
the  small  npn  transistor  of  the  741  has  a  collector  substrate  capacitance 
of  3.2  pF.  The  lateral  pnp  transistor  has  a  base  substrate  capacitance 
of  5.1  pF.  Taking  these  values  as  a  starting  point,  a  nonlinear  macro¬ 
model  for  the  741  which  includes  substrate  capacitances,  is  shown  in 
Figure  2-8.  The  true  value  of  the  individual  device  substrate  capaci¬ 
tance  in  situ  on  the  IC  is  uncertain  due  to  the  complexity  of  the  IC 
structure.  Fang  has  found  that  for  the  purpose  of  numerical  prediction, 
the  exact  value  of  substrate  capacitance  is  unnecessary.  A  typical 
capacitance  value  can  be  used  to  obtain  good  EMI  predictions.  The 
same  situation  also  is  true  for  the  nonlinear  macromodel.  Consequently, 
one. substrate  capacitance  value  is  used  at  both  the  collector  and 
the  emitter  for  both  BJTs  in  the  nonlinear  macromodel.  It  will  be 
shown  later  that  a  structural  justification  for  the  substrate  capaci¬ 
tance  values  (which  cannot  be  given)  is  not  essential.  What  is  essential 
Is  a  methodology  for  accounting  for  the  effects  of  substrate  capacitance 
upon  the  Internal  RF  interaction  which  produces  the  intermodulation. 

At  the  present  time  no  procedure  exists  by  which  the  values  for  the  para¬ 
sitic  capacitors  C$ub  shown  in  Figure  2-8  can  be  determined  independ¬ 
ently.  Instead  a  parametric  fitting  procedure  is  used.  Values  such  as 
0,  2  pF,  or  4  pF  are  assigned  to  the  four  capacitors  C$ub  and  the  NCAP 


calculations  are  made. 


The  results  of  the  NCAP  predictions  for  the  741  buffer  amplifier 
second  order  transfer  function  using  substrate  capacitance  values  of 
0  pF,  2  pF  and  4  pF  are  shown  in  Figure  2-9.  It  is  noted  that  the 
inclusion  of  substrate  capacitance  in  the  nonlinear  macromodel  improves 
the  numerical  accuracy  of  prediction  and  also  gives  better  qualitative 
agreement  with  the  experimental  results  over  a  wider  frequency  range. 

A  substrate  capacitance  of  4  pF  produces  the  best  fit  between  NCAP 
predictions  and  experimental  results.  Furthermore,  NCAP  predictions 
made  using  a  Op  Amp  macromodel  with  C^  =  4  pF  compare  favorably 
with  Fang's  NCAP  predictions  made  using  the  Op  Amp  full  model  with 
Csub  *  2  PF  or  Wooley's  C$ub  values.  See  Figure  8  *n  Appendix  A. 

When  the  macromodel  substrate  capacitance  value  is  increased  above 
4  pF,  the  second  order  transfer  function  values  predicted  change  less 
rapidly  and  tend  to  exceed  the  experimental  values. 

The  computer  time  for  an  NCAP  analysis  is  greatly  reduced  by  the 
use  of  the  nonlinear  macromodel.  On  the  Honeywell  6180  computer,  an 
NCAP  five-frequency  sweep  analysis  of  the  741  Buffer  Amplifier  second 
order  transfer  function  based  on  the  nonlinear  macromodel  requires  34 
seconds  of  CPU  time.  In  addition  30  seconds  for  peripheral  transferring 
of  data  Is  required.  Not  only  does  the  nonlinear  macromodel  reduce 
computer  time  and  costs,  it  is  also  much  easier  for  the  EMC  engineer 
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Figure  2-9.  Effect  of  the  BJT  Substrate  Capacitance 
Value  on  NGAP  Prediction  of  the  Second 
Order  Transfer  Function  Based  on  the 
Nonlinear  Macromodel  of  the  ?41  Buffer 
Amplifier. 


CHAPTER  THREE 


MACROMODELING  PROCEDURE  AND  APPLICATION 

In  the  previous  chapter  an  Op  Amp  nonlinear  small -signal  macro¬ 
model  based  upon  the  Op  Amp  nonlinear  large-signal  macromodel  of  Boyle 
et  al..^2^  was  used  to  predict  EMI  in  a  741  Op  Amp.  One  obvious  ques¬ 
tion  is  can  the  Op  Amp  small -signal  macromodel  be  used  for  other  Op 
Amp  types.  If  the  answer  is  yes,  a. second  question  is:  What  adjustments 
are  needed  for  other  Op  Amp  types?  It  is  believed  that  the  Op  Amp  non¬ 
linear  small-signal  macromodel  shown  in  Figure  2-8  is  global  in  the 
same  sense  as  the  BJT  nonlinear  T-model  shown  in  Figure  1-2  is  global. 
The  Op  Amp  nonlinear  small -signal  macromodel  can  be  used  for  any  inte¬ 
grated  Op  Amp  (including  JFET-bipolar,  MOSFET-bi polar,  and  MOSFET  types) 
when  three  types  of  adjustments  are  made: 

1.  The  Input  transistors  and  Q2  in  Figure  2-8  must  cor¬ 
respond  to  the  input  transistor  types  in  the  actual  Op 
Amp.  Appropriate  NCAP  transistor  model  parameters  for 
the  transistors  Qj  and  Q2  must  be  used. 

2.  Appropriate  values  for  the  resistors,  capacitors  and 
transconductances  in  the  original  Boyle  Op  Amp  macromodel 
must  be  determined  using  the  techniques  described  by  Boyle 
et  al.^2^  for  bipolar  Op  Amps  and  by  Krajewska  et  al.^®^ 
for  FET-bi polar  Op  Amps. 


3.  Values  for  the  four  substrate  capacitances  C$ubl ,  Csub2, 

^sub3*  an<*  ^sub4  s*10wn  fri9ure  2-8  must  be  determined 
using  the  procedures  described  in  this  chapter. 

In  this  chapteY*  the  basic  macromodeling  procedures  will  be  outlined 
in  more  detail  in  Section  3.1.  In  Sections  3.2  and  3.3  the  procedures 
will  be  applied  to  a  new  Op  Amp-the  LM10.  The  LM10  Op  Amp  is  a  bipolar 
Op  Amp,  but  one  that  is  very  different  from  the  741  Op  Amp  discussed  in 
Chapter  2.  In  Section  3.4  NCAP  EMI  predictions  and  experimental  results 
will  be  compared  in  order  to  show  how  to  determine  appropriate  values 
for  the  substrate  capacitances  shown  in  Figure  2-8. 

3.1  The  Macromodeling  Procedure 

In  this  section  the  basic  macromodeling  procedure  for  Op  Amps  is 
described.  The  specific  case  for  BJT  Input  transistors  is  considered. 

A  similar  procedure  can  be  used  for  JFET  or  MOSFET  input  transistors. 

(1)  A  BJT  IC  Op  Amp  is  represented  by  the  model  of  Figure  2-8.  Note 
that  substrate  capacitances  are  included. 

(2)  The  NCAP  BJT  model  parameters  for  the  input  BJT  pair  are  taken 
from  sources  such  as  the  work  of  Fang.^*^  when  device  fabrica¬ 
tion  data  are  available,  it  may  be  possible  to  extract  the  model 
parameters  using  computer-aid  analysis  techniques. 

(3)  The  remaining  Boyle  macromodel  parameters  are  extracted  from  the 
Op  Amp  performance  specifications  using  procedures  described 

by  Boyle  et  al.^-25^ 


(4)  The  IC  Is  tested  in  the  experimental  system  of  Figure  1-3  in  the 
unity  gain  buffer  amplifier  configuration.  The  data  are  con¬ 
verted  into  values  for  the  second  order  transfer  function 
Hgjfg.-fi )  using  Equation  (1-14).  The  second  order  transfer 
function  is  plotted  vs  frequency  as  shown  in  Figure  1-6.  The 
second  order  nonlinear  transfer  function  ^(fg.-f-j)  represents 
well  the  experimental  manifestation  of  EMI  effects  related  to 
the  demodulation  of  AM  modulated  RF  signals  in  the  Op  Amp. 

(5)  A  NCAP  computer  simulation  is  performed  to  obtain  predicted 
values  for  H2(f2,-f.j).^56l  The  substrate  capacitance  values 
are  then  adjusted  to  obtain  the  best  fit  between  predicted 
and  experimental  results. 

The  macromodeling  procedure  just  described  will  be  applied  to  the 
LM10  Op  Amp. 

3.2  Macromodeling  of  the  LM10  Bipolar  IC  Op  Amp 

r57i 

The  LM10  Op  AmpL  J  has  on  one  monolithic  chip  an  Op  Amp,  a  . 
precision  voltage  reference,  and  an  adjustable  reference  buffer 
amplifier.’  The  pin  configuration  and  functional  diagram  is  shown  in 
Figure  3-1.  The  Independent  Op  Amp  portion  has  pin  to  pin  equivalency 
with  the  741  Op  Amp.  Upon  examining  the  National  LM10  data  sheets 
(Appendix  B),^]  it  is  noted  that  the  LM10  Op  Amp  performance  is 
similar  to  that  of  the  741  Op  Amp.  However,  the  similarity  stops 
here.  The  LM10  Op  Amp  has  been  designed  so  that  it  can  be  operated 
from  a  single  supply  voltage’ as  low  as  1.1  V.  The  LM10  Op  Amp  IC 
schematic  circuit  diagram  is  shown  in  Figure  3-2. 
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Figure  3-1 • 


LMLO  Pin  Configuration  and  Functional 

_4  65% 

Diagram. 


The  Op  Amp  input  section  of  the  LMIO  is  different  from  the  741  Op 
Amp  in  major  design  respects.  The  essential  details  of  the  LM10  Op 
Amp  input  section  are  shown  in  Figure  3-3.  gf  particular  signific¬ 
ance  is  that  the  LM10  input  transistors  are  pnp  devices.  Thus,  the  LM10 
Op  Amp  constitutes,  to  a  considerable  degree,  a  new  test  on  the  effec¬ 
tiveness  of  the  nonlinear  macromodel  for  EMI  analysis.  The  macromodel 
for  the  LMIO  is  shown  in  Figure  3-4.  The  macromodel  equations  and  para¬ 
meter  value  extraction  methods  are  presented  in  the  next  section. 

3.3  Calculation  of  the  Macromodel  Parameters 


The  LMIO  macromodel  parameter  model  parameters  refer  to  the  equiv¬ 
alent  circuit  of  Figure  3-4.  In  this  section  values  for  the  para¬ 
meters  are  determined  using  the  procedures  given  by  Boyle  et  al.^29^ 
and  data  from  National  Semiconductor  data  sheets  for  the  LMIO  which  are 
included  as  Appendix  B.  For  a  complete  understanding  of  the  procedures 
used,  a  reading  of  Boyle's  paper  is  essential.  For  the  input  pnp  BJTs, 
the  dc  collector  current  I£  is  1.25  yA  according  to  Wildar  et  alJ59^ 
The  transconductance  g^  of  the  input  differential  stage  is  given  by^] 


.  _  qIC  _  1.25  x  ID'6.  „  ..  ,„-5  .... 

irnr'  2  x  0.0258  “  2*42  x  10  A/v  • 

The  compensation  capacitance  C ^  is  given  by 

*  “  ^  -  2-42  x  10  5  =  6.4  x  10"11  F 


(3-1) 


Z*f0dB  2 tt  x  60  x  10' 


(3-2) 
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The  frequency  fQdB  at  which  the  LM10  Op  Amp  open-loop  gain  is  0  dB 
Is  60  kHz.  The  resistor  »  R5  +  R6  *  R8  *  1G4  kQ  is  taken  from 
Figure  3-2.  The  conductance  Ga  is  the  reciprocal  of  R£. 


1 

5 


1  -j  *  9.6  x  IQ"6  A/V 


(3-3) 


a  ^  104  x  10^ 

Next  the  output  resistances  Rg-j  and  Rgg  are  calculated.  The  data 
sheets  (Appendix  B)  show  an  output  Impedance  varying  with  frequency  as 
a  function  of  load  current  over  four  orders  of  magnitude  for  closed-loop 
gains  of  1  and  100.  It  is  possible  to  use  these  curves  to  determine  the 
Op  Amp  open-loop  output  resistance.  The  output  resistance  R.  t  =  + 

Rq2  can  also  be  determined  from  the  open-loop  dc  voltage  gaiVi  (AvD)  vs 
load  resistance. (R^)  graph  shown  in  Appendix  B  using  the  relationship 

Avd  *  Ayp  (R^  *  *)  x  R^/ (Rl  +  Rjnt)  (3-4) 


where  AvQ  (R^  *  «)  is  the  open-circuit  value  for  AyD.  Choosing  points 
R^  ■  2  kfl  where  AvQ  *  110  dB  (3.16  x  105)  and  R^  =  10  kn  where  Ay[)  = 
114  dB  (5.01  x  1 05)  and  using  Equation  (3-4),  the  value  determined  for 
R1nt  ^13  The  macromodel  dc  output  resistance  parameter  RQ1  is 
arbitrarily  set  at  0  a.  The  macromodel  ac  output  resistance  parameter 
Rq2  Is  set  equal  to  R^^  *  1713  n. 

Hooley's^-24^  value  for  the  dc  current  gain  Bp  for  a  lateral  pnp  is 
0p  *  75  which  is  close  to  the  value  Bp  =  100  given  by  Widlar.^3^  The 
Hooley  value  is  used.  The  value  for  the  Early  Voltage  given  by  Fang 


In  Table  4-2  of  Reference  18  for  the  lateral  pnp  transistors  Q3  and 
Is  VA  *  55  V.  This  value  is  used.  Using  I£E  =*  2  Iq,  the  resistance 
R£  is  given  by 

Re  *  VA/IEE  =  55/ (2. 5  x  10"6)  *  2.2  x  107  a  (3-5) 

The  capacitance  is  calculated  from  the  frequency  response  curves  for 
the  gain  and  the  phase.  From  the  frequency  response  curves  in  the 
specifications  (Appendix  B),  the  phase  margain  <|>m  of  the  open-loop  re¬ 
sponse  is  obtained  at  the  frequency  where  the  gain  becomes.  0  dB.  The 

*  value  obtained  is 
m 

a  «  180°  -  120°  =  60° 
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The  excess  phase  shift  at  f^g  as  defined  by  Boyle^29^  is  then, 

A+  *  90°  -  60°  =  30° 

The  excess  phase  shift  is  related  to  C-j  by 

C1  *  (C2/2)tanA<|>  -  (64/2)tan30°  =  18.5  pF  (3-6) 

The  value  for  hpE  is  calculated  by  taking 

hFBnax  *  100’  and 

hFE  “  hFEmax/{1  +  alog2^IC/ICmax^}  *3-7^ 


100/{1  +  0.8071ogZ(l .25/110)}  *  24.7 


where  the  value  a  *  0.807  Is  the  value  given  by  Fang  for  the  lateral  pnp 
transistors  Q3  and  Q4  In  Table  4-3  of  Reference  18.  The  emitter  re¬ 
sistance  Is  selected  to  make  the  gain  of  the  input  differential 
amplifier  stage  ^/v^  *  1.^^  Setting  vd/v.jn  *  1,  Boyle  obtained  the 
expression 

Re  -  [Rc  '  l/gJ[hFE/(hFE  ♦  U1  (3-8) 

Using  Equation  (3-8)  and  the  values  R£  *  104  kfl,  gm  *  qIc/kT  *  48.4  pA/V, 
and  hFE  *  24.7,  the  value  Rfi  »  80  kn  Is  obtained. 

The  transconductance  is  calculated  using  the  expression  for  the 
low  frequency  open-loop  gain  ayD  given  by  Boyle^-29^ 

*VD  "  (VvJ(GaR2><GbR02>  <3-9) 

The  value  R2  *  100  kn  used  by  Boyle^29^  for  the  741  Op  Amp  is  also  used 
for .the  LM10.  Using  Equation  (3-9)  with  a  typical  value  avD  =  300,000, 
v^v.^  ■  1,  Rg  *  100  kn,  Gfl  »  9.6  pA/V,  and  Rq2  *  1713  n,  the  value 
■  182.4  A/V  is  calculated. 

The  capacitance  is  related  to  the  negative  going  slew  rate  by 
the  formula 

SJ  -  (2Ic)/(C2  +  CE)  (3-10) 

Using  the  graphs  shown  In  Appendix  B  for  the  comparator  response  time 
for  various  input  drives,  it  is  observed  that  the  output  voltage  changes 
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by  5  V  In  0.2  ms  which  corresponds  to  a  slew  rate  =  5  V/0.2  ms  = 

2.5  x  10*  V/s.  Using  Equation  (3-10)  with  =  2.5  x  10*  V/s  and 
Cg  *  64  pF,  the  value  36  pF  is  obtained  for  C^. 

The  lateral  pnp  transistor  NCAP  parameter  values  determined  using 
Wooley's  data^24^  are  used  for  the  transistors  and  Q2  in  the  macro¬ 
model.  These  values  are  those  for  transistor  Q3#  in  Table  1-3. 

3.4  NCAP  Predictions  for  the  Second  Order  Transfer  Function  H2(f2,-f1) 

Based  on  the  Nonlinear  Macromodel 

The  NCAP  coding  diagram,  using  the  parameters  given  in  the  previous 
section,  is  shown  in  Figure  3-5.  The  input  code  for  an  NCAP  calculation 
of  the  second  order  nonlinear  transfer  function  H2(f2,-f-| )  is  listed  in 
Table  3-1.  The  NCAP  predicted  results  presented  in  Figure  3-6  do  not 
agree  well  with  the  experimental  results  because  substrate  capacitance 
effects  are  not  accounted  for. 

The  actual  location  for  the  lateral  pnp  transistor  parasitic  sub¬ 
strate  capacitance  is  from  the  base  terminal  to  ground  (substrate)  as 
shown  In  Figure  3-7. ^2*-l  The  lateral  pnp  transistor  in  the  Op  Amp  in¬ 
vestigated  by  Wooley  has  a  base  area  approximately  4.5  mil  x  3  mil  and 
a  substrate  capacitance  of  1  pF.  A  measurement  of  the  LM10  die  area  for 
the  base  of  Q2  yields  a  value  approximately  3  mil  x  6  mil.  This  area  is 
nearly  equal  to  that  of  the  lateral  pnp  transistors  Q3  and  Q4  listed  in 
Table  1-3.  Notwithstanding  the  difference  in  the  actual  physical  place- 


TABLE  3-1 

NCAP  INPUT  CODE  FOR  THE  MACROMODEL  LMIO  VOLTAGE  FOLLOWER 


**  MACROMODEL  FOR  LMIO  OP  AMP 

*♦  ANALYSED  BY  NCAP  FOR  H2 

?*  FOR  STANDARD  SUPPLY 

**  VOLTAGE  OPERATION 

**  IN  BUFFER  AMPLIFIER 

♦*  UNITY  GAIN  CONFIGURATION 

♦*  SUBMITTED  AT  SUNY  BUFFALO 

**  CCS«0  .  CES*0  .  CBS=0 

**  OUTPUT  NODE  IS  15  PIN  6  IS  NODE  14 

*  START  CIRCUIT 

*  GENERATOR 
NODE  1  0 

FR  1  1E7  8E7  5  LIN 
AMP  1  0.0 

FR  2  -0.*99996E7  -7.99996E7  5  LIN 
AMP  1  0.0 
IMP  50  0 

♦LINEAR  COMPONENTS 
R  2  1  0.01 
R  2  3  560 
R  3  8  4E3 
R  2  0  50 
R  0  6  104E3 
R7  If  80E3 
R  0  10  104E3 
R  11  16  80E3 
R  16  0  22E6  ' 

R  12  0  1.0E5 
R  13  0  1713 
R  13  14  0.1 
R  14  15  4.7E4 
R  15  0  4.489E4 
R  14  17-  620 
R  17  4  4E3 
C  2  0  1.1E-10 


C  6  10  18.5E-12 
C  16  0  36E-12 
C  12  13  64E-12 
C  15  0  2.106E-9 

**  SUBSTRATE  CAPACITANCES  CES  =  0  PF, 
**  CBS  =  0  PF,  CCS  =  0  PF 
C  4  0  5.1E-15+ 

C  8  0  5.1E-15 
C  6  0  4E-15 
C  10  0  4E-15 
C  7  0  0.5E-15 
C  11  0  0.5E-15 
♦TRANSISTOR 
NODE  4 

1.5  12.6  95  0.333  1.25E-6 
0.11E-3  0.808  100  0.834E-12 
1.0  0.1E-12  1.06E-6  500 
0.5E6  0.1E-12  0.1E-12 
♦TRANSISTOR 
NODE  8 

1.5  12.6  95  0.333  1.25E-6 

0.1 IE- 3  0.807  100  0.834E-12 

1.0  0.1E-12  1.06E-6  500 

0.5E6  0.1E-12  0.1E-12 

♦LINEAR  DEPENDENT  SOURCE 

NODE  6  10  0  12 

VC  9.6E-6  0.0 

♦LINEAR  DEPENDENT  SOURCE 

NODE  16  0*120 

VC  0.79E-11  0.0 

♦LINEAR  DEPENDENT  SOURCE 

NODE  12  0  0  13 

VC  182.4  0.0 

♦END  CIRCUIT 

♦END 


t  For  ease  of  adding  and  eliminating  substrate  capacitances  in  NCAP 
runs,  a  0  pF  capacitance  is  represented  by  a  reduction  of  three 
orders  of  magnitude. 
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merit  of  the  substrate  capacitance  from  base  to  ground  as  shown  in 
Figure  3-7,  the  second  order  transfer  function  fitting  procedure 
outlined  in  Section  3.1  requires  that  an  emitter  to  ground  substrate 
capacitance  be  used.  That  this  substrate  capacitance  location  does  not 
correspond  to  a  physical  location  in  the  LM10  Op  Amp  implies  that  the 
emitter  substrate  capacitance  i?  an  equivalent  parameter.  The  second 
order  transfer  function  l^^t-f-])  is  very  sensitive  to  any  change  in 
this  capacitance,  as  illustrated  in  Figure  3-8.  On  the  other  hand  a 
non-zero  base-substrate  capacitance  Cg^  which  has  more  physical  signif¬ 
icance  reduces  the  accuracy  of  the  nonlinear  macromodel  predictions  as 
Illustrated  in  Figure  3-9.  This  apparent  contradiction  may  be  a  re¬ 
sult  of  simulating  the  many  distributed  substrate  parasitic  capaci¬ 
tances^58^  with  a  few  discrete  capacitors  placed  at  the  input  tran¬ 
sistor  terminals. 


The  final  choice  for  the  LM10  nonlinear  macromodel  substrate 
capacitances  obtained  by  comparing  NCAP  predicted  and  experimental 
values  for  the  second  order  transfer  function  H2(f2,-fj)  are  CgS  *  0  pF 
Ccs  *  0  pF  and  *  0.5  pF.  The  NCAP  calculated  values  for  the  first 
order  transfer  function  H^(f)  for  the  LM10  unity  gain  buffer  amplifier 
are  compared  to  experimental  results  in  Figure  3-10.  The  NCAP  predic¬ 
tions  and  experimental  results  for  the  second  order  transfer  function 
82^2’~^1 )  are  shown  in  Figure  3-11. 


As  shown  in  Figure  3-11  by  properly  adjusting  the  macromodel  sub¬ 
strate  capacitance  values  excellent  agreement  has  been  obtained  between 
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Transfer  Function  of  LM10  Buffer  Amplifier. 


NCAP  predicted  values  and  experimental  values  for  Hgtfg.-fj)  for 
frequencies  above  3  MHz.  For  frequencies  below  1  MHz  the  predicted 
values  exceed  the  measured  values  by  as  much  as  18  dB  (at  300  kHz). 
Much  of  the  difficulty  In  achieving  good  agreement  at  frequencies  less 
than  1  MHz  may  be  a  result  of  the  low  H2(f2,-f^)  values  (-35  dB) 
measured  below  1  MHz.  A  similar  difficulty  was  encountered  using  a 
full  model  for  the  yA741  Op  Amp.  (See  Figure  1-11.) 

It  Is  Interesting  to  compare  yA741  and  LM10  Op  Amp  experimental 
values  for  H2(f2,-fj).  As  shown  in  Figure  3-12  at  frequencies  above 
1  MHz  the  H2(f2,-f1)  experimental  values  for  the  LM10  Op  Amp  are  at 
least  15  dB  below  those  for  the  yA741  Op  Amp.  This  means  that  the 
LM10  Op  Amp  Is  less  susceptible  to  EMI  at 'frequencies  above  1  MHz  than 
the  yA741  Op  Amp.  The  lower  LM10  Op  Amp  susceptibility  to  EMI  above 
1  MHz  manifests  itselfs  in  the  lower  Csufa  values  used  in  the  LM10  Op 
Amp  macromodel.  Lowering  the  values  lowers  the  predicted 
H2(f2*-fi)  values. 

The  reason  for  the  lower  LM10  Op  Amp  susceptibility  to  EMI  may  be 
Its  total  compensation  capacitance  value  of  1,000  pF.^58^  This  large 
capacitance  is  achieved  by  the  use  of  emitter  isolation  junction 
capacitances  of  1  pF/mil  .  Also  the  intermediate  stage  has  a  phase 
compensation  capacitance  of  100  pF  (Cl  of  Figure  3-3)  and  an  overall 
compensation  capacitance  from  the  output  to  the  positive  input  tran¬ 
sistor  of  22  pF,  (C2  of  Figure  3-3).  By  performing  a  full -model  NCAP 
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Experimental  Values  for  the  Second  Order  Transfer  Function 
vs  Fluency  for  the  741  and  LM10  Op  Amps. 


analysis  of  the  LM10  Op  Amp,  It  Mould  be  possible  to  do  a  parameter 
sens1t4vity  analysis  to  determine  the  importance  of  the  various  capacit 
ors  in  the  LM10  Op  Amp.  This  has  not  been  attempted  because  the  goal 
of  the  present  effort  is  not  to  evaluate  Op  Amp  full  models  but  to 
evaluate  Op  Amp  macromodels.  The  ability  of  the  Op  Amp  macromodel  to 
account  for  EMI  in  the  fairly  complex  integrated  LM10  Op  Amp  without 
resorting  to  a  full  model  analysis  is  the  important  point. 


CHAPTER  FOUR 


NETWORK  THEORY  FOR  THE  NONLINEAR  MACROMODEL 

In  the  two  previous  chapters.  It  has  been  demonstrated  that  a 
nonlinear  small-signal  Op  Amp  macromodel  which  contains  only  two  BJTs 
Is  effective  and  accurate  for  EMI  analysis.  Furthermore.  It  Is  believed 
that  the  nonlinear  macromodel  shown  In  Figure  2-2  Is  global  and  can 
be  applied  to  all  bipolar  Op  Amps  of  present  technology  and  design*  An 
actual  bipolar  Op  Amp.  such  as  the  one  Illustrated  In  Figure  1-8,  Is  of 
course  made  of  many  Identifiable  (and  some  unidentifiable)  components. 

All  the  BJTs  In  the  Op  Amp  have  nonlinear  parameter  values  which  are 
comparable,  as  shown  by  the  work  of  Fang.^18^  How  then  is  it  that  the 
RFI  contributions  of  the  large  number  of  nonlinear  components  in  the 
Op  Amp  can  be  represented  by  just  two  nonlinear  transistors  in  the 
macromodel?  This  question  is  addressed  in  the  present  chapter. 

4.1  Cascade  Configuration  for  the  IC  Operational  Amplifier 

A  common  IC  Op  Amp  design  Is  the  three  stage  configuration  utilized 
T421 

In  the  uA741.L  J  This  configuration  consists  of  an  input  stage,  an 
Intermediate  gain  stage  with  compensation,  and  a  unity  gain  output  stage. 
For  small-signal  analysis  purposes,  the  IC  schematic  diagram  can  be 
simplified  to  that  of  Figure  4-1.  A  small-signal  linear  equivalent 
circuit  shown  In  Figure  4-2  Illustrates  the  three  stage  cascade  more 
clearly. 


Input  Stage 


Gain  Stage 


Output  Stage 


Figure  4-2. 


Small-Signal  Equivalent  Circuit  of  the 
IC  Op  Amp. 


For  small -signal  nonlinear  analysis,  each  section  in  Figure  4-2  is  re¬ 
placed  by  an  appropriate  small-signal  nonlinear  representation.  One 
such  representation  is  shown  in  Figure  4-3  in  which  A^f)  is  the  first 
order,  or  linear,  transfer  function,  A^f-j.fg)  Is  the  second  order 
transfer  function  and  A3(f1,f2,f3)  is  the  third  order  transfer  function. 
When  needed,  higher  order  transfer  functions  can  also  be  included.  The 
summing  node  indicates  that  all  transfer  functions  operate  on  the  input 
signal  frequency  spectrum  U(f)  to  produce  the  total  output  response 
spectrum  V(f). 

A  condensed  block  diagram  for  the  Op  Amp  three  stages  is  shown  in 

* 

Figure  4-4(a)  as  a  three-stage  cascade  .  The  three-stage  cascade  can 
be  viewed  as  two  two-stage  cascades  as  illustrated  in  Figures  4-4(b) 
and  (c).  Therefore  only  the  cascading  relationship  for  two  stages  need 
be  considered.  The  block  diagram  shown  in  Figure  4-3  can  be  substituted 
Into  each  stage  of  the  two-stage  cascade  shown  in  Figure  4-4(b)  to  ob¬ 
tain  the  block  diagram  shown  in  Figure  4-5.  The  two-stage  cascade 
shown  in  Figure  4-5  has  transfer  functions  D^f^),  D2(f1,f2),  and 
D3(fi,f2,f3)  which  can  be  written  in  terms  of  the  single-stage  transfer 
functions  using  the  cascade  relationships  given  by  Weiner  and  Spina**. 

*  This  first  order  approximation  of  the  2-port  cascading  system  as¬ 
sumes  the  loading  effect  of  the  stages  to  be  negligible. 

**  A  derivation  of  the  cascade  relationships  according  to  Weiner  and 
Spina  is  given  in  Chapter  4  of  Reference  15. 


Figure  4-4.  Diagrammatic  Representation  of  Cascaded 
Stages  Making  Up  an  Op  Amp. 


(a)  Three-Stage  Representation  of  the 
Op  Amp  Equivalent  Circuit. 

f 

(b)  Two-Stage  Cascade  for  Transfer 
Function  Derivation. 

(c)  Conversion  of  the  Three-Stage  System 
into  a  Two-Stage  System. 


Figure  4-5 •  A  Two-Stage  Cascade  Nonlinear  Transfer 
Function  Representation. 


Figure  4-6.  Block  Diagram  Illustrating  Nonlinear  Transfer 
Functions  of  Op  Amp  Input  Stage  and  Second 
Stage  Cascade  With  Feedback  in  the  Second  Stage 


The  results  are: 


D-j  (f-j )  -  B-|  (f-j  )A-j  ( f -j )  (4-1) 

D2(frf2)  -  B1(f1+f2)A2(f1  ,f2)  +  B2(fltf2)?  A^)  (4-2) 

°3  (f1,f2»f3)  *  Bl(fl rf2+f2)A3(fl*f2*f3) 

♦  2B2(flff2+f3)A1(f1)A2(f2,f3) 

♦  B3(fvf2,f3)|  A^)  (4-3)* 

Consider  the  second  order  transfer  function  of  two  cascaded  stages 
which  each  have  comparable  nonlinear! ties.  Equation  (4-2)  Implies  that 
the  contributions  of  both  the  first  and  the  second  stage  to  the  two- 
stage  cascade  second  order  transfer  function  need  to  be  considered. 
Similarity,  Equation  (4-3)  Implies  that  the  contributions  of  each  stage 
to  the  two-stage  cascade  third  order  transfer  function  need  to  be  con¬ 
sidered.  This,  observation  Is  supported  by  the  experiments  of  Narayanan^- 
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and  Weiner  and  Sp1nak  J  Involving  discrete  BJT  cascades.  However,  the 
feedback  internal  to  the  IC  Op  Amp  strongly  affects  the  two-stage 
cascade  second  order  transfer  function  02(f^,f2)  when  the  sum  frequency 
fj  +  f2  Is  low  enough.  This  point  will  be  discussed  in  the  next  section 


*  See  Reference  15,  Chapter  4  for  the  procedure  of  averaging  which 
makes  up  the  overbar  term  In  Equation  4-3. 


4.2  The  Effect  of  Internal  Feedback  on  A  Cascade  Nonlinear  Transfer 
Function 


The  Internal  feedback  capacitor  C2  shown  in  the  small-signal  equiv¬ 
alent  circuit  of  Figure  4-2  has  a  major  effect  at  RF  frequencies.  To 
emphasize  the  Important  role  played  by  the  Internal  feedback  capacitor 
C2  the  Input-stage-gain-stage  cascade  for  the  Op  Amp  is  drawn  as  shown 
In  the  block  diagram  of  Figure  4-6,  where  G  represents  the  input  stage 
transfer  functions,  H  represents  the  gain  stage  transfer  functions,  and 
B  represents  the  feedback  of  the  internal  compensation  capacitor  C2. 
Employing  Narayanan's  algorithm  (see  Equations  8  and  11  of  Reference  8), 
the  first  order  transfer  function  F^(f)  of  the  second  stage  with  the 
linear  feedback  e  is:  - 


F,(f) 


H,(f) 


(4-4) 


and  the  second  order  transfer  function  F2(fj,f2)  of  the  second  stage 


with  the  feedback  Is: 


H2(frf2)  4  [1  -s(f i ) F-j  (f .)] 

F2(frf2}  "  1  +  6(ii1+f2)H1(f1+f2)  ' 


(4-5) 


Substituting  Equation  (4-4)  into  Equation  (4-5)  leads  to  the  result 

„  /Jf  „  ,  H2(frf2) 

F2(frf2}  a  Ll+H1(f1)Blf1)ILl+H1(f2)B(f2)]Cl+H1(f1+f2)B(f1+f2)] 

(4-6) 

Note  that  Equation  (4-6)  for  the  second  order  transfer  function  F2^fl’f2^ 
has  been  written  in  terms  of  the  first  order  transfer  function  H.(f). 


The  gain  stage  of  the  yA741  Op  Amp  which  Is  now  considered  as  a  specific 
case  has  the  small-signal  equivalent  circuit  shown  In  Figure  4-7.  The 
feedback  factor  3  Is  given  by 

6  -  1/Cl  -  j/2*fRC2]. 

For  R  *  100  kfl  and  »  30  pF,  Equation  (4-7)  can  be  written  as 

3  -  1/Cl  -  j*05/f(MHz}] 

since  the  cutoff  frequency  l/2irRC,,  »  .05  MHz.  The  feedback  factor  3 
has  a  magnitude  0.008  at  400  Hz,  0.7  at  0.050  MHz  and  1.0  at  frequencies 
greater  than  0.2  MHz.  Typical  values  for  the  open-loop  second-stage 
first  order  voltage  gain  H^(f)  are  estimated  to  be  30,000  at  400  Hz, 

250  at  .050  MHz  and  less  than  0.1  at  9  MHz.  The  denominator  of  Equation 
(4-6)  varies  over  a  wide  range  of  values.  Let  us  consider  the  specific 
case  corresponding  to  the  demodulation  of  AM  modulated  RF  signals  by 
the  second  stage.  The  frequency  f^  +  f2  Is  kept  constant  at  400  Hz. 

The  frequency  f2  Is  varied  over  the  range  .05  to  100  MHz.  Ts:e  frequency 
f-j  Is  negative  and  must  track  the  frequency  f2  so  that  the  condition 
fj  +  f2  »  400  Hz  Is  satisfied.  The  factor  [1  +  H-j (f1+f2)3(f-|+f2)] 
has  a  value  given  by 

[1  ♦  H1(f1+f2)3(f1+f2)]  »  [1  +  (30 ,000) (0. 008)]  =  240  (4-9) 

which  does  not  vary  as  the  frequency  f2  Is  varied  because  fj  +  f2  Is* 
kept  constant  at  400  Hz.  At  a  low  RF  frequency  such  as  f2  *  .050  MHz, 


(4-7) 


(4-8) 
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the  factors  [1  +  (f ^ )e(f -j )3  and  [1  +  H-j  (f2)e(f2)]  are  of  the  order 

[1  +  (250) (0.7)3  *  200  and  the  denominator  of  Equation  (4-6)  is  of  the 
order  of  (200) (200) (250)  *  TO7.  Clearly  at  an  RF  frequency  f2  ■  0.050 
HHz  the  second-stage  second  order  transfer  function  F2(f^,f2)  is  reduced 
essentially  to  zero  by  the  feedback  Internal  to  that  stage.  For  an  RF 
frequency  f2  above  1  MHz,  the  factors  [1  +  H-j (f-j )e(f-j )]  and 
[1  +  Hi (f2)e(f 2)3  have  a  value  near  unity  and  the  denominator  of 
Equation  (4-6)  has  a  constant  value  240  given  by  Equation  (4-9).  For 
an  RF  frequency  f2  >  1  MHz,  the  open-loop  second-stage  second  order 
transfer  function  H2  (f^,f2)  Is  reduced  by  a  factor  240  by  the  feedback 
Internal  to  that  stage.  Thus  the  contribution  of  the  second-stage  to 
the  Op  Amp  total  second-order  transfer  function  is  suppressed  signific¬ 
antly  by  the  internal  feedback  capacitor  C2  of  that  stage.  This  suggests 
(but  does  not  prove  absolutely)  that  the  nonlinear  effects  associated 
with  the  second  stage  can  be  neglected.  It  is  believed  that  similar 
analysis  would  show  that  the'thlrd  order  transfer  function  F3(f.j 
of  the  second  stage  Is  also  reduced  by  the  internal  feedback  capacitor 
C2.  Thus  It  seems  reasonable  to  simplify  the  Op  Amp  first-stage  and 
second-stage  cascade  representation  from  that  shown  in  Figures  4-5  and 
4-6  to  that  shown  In  Figure  4-8.  In  Figure  4-8  the  second-stage  second 
and  third  order  transfer  functions  are  omitted. 

Now  the  Op  Amp  output  stage  needs  to  be  considered.  As  shown  in 
Figure  4-1,  the  output  stage  is  a  unity  gain  stage  consisting  of  com- 


84 


plementary  emitter  followers,  knitter  followers  may  be  viewed  as  feed¬ 
back  amplifiers  having  Internal  unity  feedback  and  a  high  loop  gain 
which  causes  the  closed-loop  gain  to  be  one.  The  analysis  and  arguments 
presented  for  the  second  stage  with  Its  Internal  feedback  should  apply 
completely  to  the  output  stage  with  Its  Internal  feedback.  Thus  the 
output  stage  second  and  third  order  transfer  functions  can  be  Ignored, 
and  the  output  stage  can  be  represented  by  a  linear  gain  stage. 

A  numerical  evaluation  of  the  second-stage  and  third-stage  contri¬ 
butions  to  the  Op  Amp  second  order  transfer  function  could  be  carried 
out  by  using  a  complete  Op  Amp  model  NCAP  analysis.  The  first  and 
second  order  transfer  function  for  the  input  stage,  the  gain  stage  with 
Cg  *  0  and  C2  *  30  pF.  and  output  stage  could  be  determined  separately. 
Then  the  cascade  relationship  expressed  by  Equation  (4-2)  could  be 
used  to  determine  how  Important  the  second-stage  contribution  is  rela¬ 
tive  to  the  first-stage  contribution.  A  similar  expression  to  Equation 
(4-2)  could  also  be  used  to  compare  the  third-stage  contribution 
relative  to  the  combined  Input-stage  gain-stage  contribution.  To  date 
such  a  numerical  evaluation  has  not  been  performed.  The  estimates  made 
using  analytical  expressions  provide  the  sole  basis  for  claiming  that 
the  Op  Amp  overall  nonlinearity  can  be  represented  by  the  nonlinearity 
of  the  Input  stage  alone  and  that  the  gain  stage  and  the  output  stages 
can  be  represented  by  linear  transfer  functions..  It  now  remains  to  be 
shown  that  the  input  stage  nonlinearity  can  be  adequately  modeled  by 
two  transistors. 


4.3  Nonlinear  Modeling  of  the  Input  Stage. 

The  Input  stage  of  the  741  is  repeated  In  Figure  4-9.  It  Is  seen 
that  some  12  BJTs  are  Involved  each  having,  nonl Inearl  ties  of  comparable 
value.  The  fact  that  the  first  stage  can  be  effectively  modeled  with  a 
differential  pair  of  transistors  In  common-emitter  (CE)  configuration 
will  now  be  discussed. 

The  operation  of  an  IC  Op  Amp  such  as  the  pA741  has  been  treated 
In  a  number  of  references. ,52^  The  transistors  Q8  through  Q12 
are  used  for  dc  biasing.  In  a  small-signal  equivalent  circuit,  they 
can  be  represented  by  linear  resistors.  The  transistors  Q1  through  Q7 
are  Involved  In  active  small-signal  processing,  and  therefore  contribute 
to  the  nonlinear  transfer  functions.  The  small-sl'gnal  equivalent  circuit 
Is  a  differential  amplifier.  Each  branch  of  the  differential  pair  Is  a 
composite  BJT  amplifier.  On  the  positive  Input  side  the  transistors 
Q1  and  Q3  form  a  common-collector  (CC)  to  common-base  (CB)  cascode. 

The  transistors  Q5  and  Q7  form  the  active  load  for  this  cascode 
amplifier.  Similarly  the  transistors  Q2  and  Q4  form  a  CC  to  CB  cascode 
with  transistor  Q6  the  active  load.  For  the  purpose  of  nonlinear  small 
signal  modeling,  only  the  cascode  amplifiers  need  be  considered. 

Modeling  the  BJTs  with  the  nonlinear  T-model  of  Figure  1-2,  the  cascode 
amplifier  formed  by  transistors  Q1  and  Q3  or  by  the  transistors  Q2  and  Q4 


Figure  4-9.  Full  Model  of  the  ?4l  Input  Stage. 
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Is  represented  by  the  small-signal  nonlinear  equivalent  circuit  of 
Figure  4-10.  An  examination  of  this  circuit  shows  that  a  three- 
terminal  representation  of  the  cascode  Is  possible.  With  terminal  C3 
being  the  composite  collector.  Cl  the  emitter  and  B1  the  base,  a 
common-emitter  stage  results  for  each  cascode  pair.  Thus  the  Input 
stage  of  the  IC  Op  Amp  can  be  represented  by  a  differential  amplifier 
formed  by  two  CE  stages.  Note  that  for  the  cascode  pair  shown  in  Figure 
4-10,  Klrchhoff's  Law  applied  to  node  El (E3)  yields 

kl(VV3>  +  rel(»2-V3)  -  k3(V4-V3)  *  Ye3(V4-V3)  (4-10) 

For  the  purpose  of  converting  the  CC  -  CB  cascode  shown  in  Figure 
4-10  to  a  single  CE  amplifier.  It  Is  not  necessary  to  solve  Equation 
(4-10).  Instead,  the  cascading  relationship  of  Equation  (4-2)  can  be 
applied.  It  may  be  noted  that  the  first  order  gain  of  the  CC  stage 
Is  unity.  The  first  order  gain  of  the  CB  stage  Is  much  greater  than 
unity.  It  also  was  reported  that  the  CB  connected  transistor  has  a 
second  order  transfer  function  considerably  lower  than  that  of  the 
CC  connected  BJT.^7*12^  Therefore,  In  applying  Equation  (4-2)  to 
Figure  4-10,  we  have  B^  high  and  A^  low  and  B^fj.fg)  <  A2(f-j,f2). 
Equation  4-2  then  can  be  written  as 

D2(fi,f2)  *  B-j (f-|+f2)  A2(frf2)  (4-11) 

which  may  be  viewed  as  second  order  transfer  function  of  a  single  stage 
amplifier. 


In  summary.  It  has  been  argued  effectively  In  this  chapter  that 
the  three-stage  amplifier  design  that  makes  up  the  741  type  of  1C  Op 
Amp  has  its  overall  nonlinear  transfer  function  dominated  by  the  input 
stage.  The  gain  stage  and  output  stage  that  follow  can  be  represented 
by  a  linear  transfer  function.  In  the  input  stage,  the  nonlinear 
transfer  function  of  CC-CB  cascode  with  active  load  can  be  represented 
by  a  single  CE  stage.  These  arguments  provide  considerable  justifica¬ 
tion  that  the  nonlinear  macromodel  of  Figure  2-2  can  be  used  for  Op 
Amp  second  order  transfer  function  calculations. 


CHAPTER  FIVE 


CONCLUSIONS  AND  RECOMMENDATIONS 

5.1  Conclusions 

This  dissertation  describes  a  small-signal  nonlinear  macromodel 
for  the  bipolar  IC  Op  Amp  which  can  be  used  successfully  for  RFI 
analysis.  Two  BJT's  are  adequate  to  characterize  the  complete  IC. 

The  small -signal  nonlinear  macromodel  can  be  obtained  directly  from 
the  large-signal  macromodel  developed  by  Boyle  et  al.^29^  The 
effectiveness  of  the  small -signal  nonlinear  macromodel  Is  verified  by 
the  agreement  between  NCAP  predictions  and  experimental  results  for 
the  second  order  transfer  function  of  the  741  Op  Amp  unity  gain  buffer 
amplifier.  Good  agreement  Is  also  obtained  between  NCAP  calculations 
based  on  the  macromodel  and  the  full  model.  By  using  the  Op  Amp  macro¬ 
model*  a  saving  of  nearly  an  order  of  magnitude  In  computer  costs  Is 
achieved.  Not  only  does  the  macromodel  conserve  computer  resources. 

It  makes  possible  for  the  EMC  engineer  to  analyse  electronic  systems 
Involving  many  Op  Amps. 

Another  Important  aspect  of  the  small -signal  nonlinear  macromodel 
Is  that  it  Is  a  global  model  which  can  be  used  for  all  Op  Amp  types. 

The  two  macromodel  transistors  should  be  the  same  type  as  those  at 
the  Input  of  the  IC  Op  Amp  being  modeled.  The  need  for  laboratory 


characterization  of  the  Individual  transistors  in  each  new  bipolar  Op 
Amp  appears  not  to  be  required.  A  combination  of  manufacturer's 
specifications  and  Fang's^18^  IC  BJT  results  can  be  used.  This  was 
demonstrated  by  the  successful  application  of  the  nonlinear  macromodel 
to  the  new  LM10  IC  Op  Amp  for  RFI  analysis  by  NCAP,  without  resorting 
to  full  model  analysis,  nor  experimental  device  characterization. 

However,  the  accuracy  of  the  nonlinear  macromodel  relies  on  the 
Inclusion  of  capacitors  to  represent  the  effects  of  substrate  parasitic 
capacitances.  The  proper  selection  of  substrate  capacitors  does  not 
necessarily  follow  the  structural  IC  parasitic  capacitances.  The  effects 
caused  by  these  capacitances  are  simulated  by  four  capacitors  in  the 
macromodel.  The  best  values  for  these  capacitors  were  found  by  compar¬ 
ing  NCAP  predicted  values  with'  experimental  values  for  a  second  order 

transfer  function  of  the  Op  Amp  In  an  unity  gain  buffer  amplifier.  Thus 
m  % 
a  parametric  fitting  Is  Involved  when  a  new  Op  Amp  Is  modeled. 

9 

That  the  combined  nonl inearl ties  of  tens  of  BJTs  In  a  linear  IC 
Op  Amp  can  be  represented  quite  well  by  a  pair  of  BJTs  In  a  macromodel 
Is  explained  by  the  cascading  theory  of  nonlinear  transfer  functions. 

The  essential  point  appears  to  be  that  the  second  stage  (gain-stage) 
contribution  to  the  Op  Amp 'second  order  transfer  function  is  reduced 
greatly  by  the  Internal  feedback  capacitor  of  the  second  stage.  Also 
It  was  argued  that  the  Input  stage  which  consists  of  two  cascode  stages 
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In  a  differential  pair  arrangement  can  be  represented  by  two  common 
emitter  stages  In  a  differential  pair  configuration  used  in  the  Op  Amp 
macromodel.  The  validity  of  the  nonlinear  macromodel  Is  therefore 
supported  both  by  experimental  verification  and  by  network  theory. 


5.2  Recommendations 


The  following  topics  are  suggested  for  future  el  rts. 

(1)  Nonlinear  macromodels  for  FET-Bipolar  IC  Op  Amp 
Large-signal  macromodels  for  JFET-Bipolar  and  M0~.  ..f-Bi polar  Op 
Amps  have  been  achieved  recently. A  small -signal  nonlinear 
macromodel  for  the  JFET-Bipolar  Op  Amp  can  be  developed  using 
the  NCAP  JFET  small -signal  nonlinear  model  now  available.  I-61-! 

A  small-signal  r.anl Inear  macromodel  for  a  MOSFET-Bi polar  Op  Amp 
can  be  developed  with  a  four-terminal  MOS  transistor  nonlinear 
model  which  Is  expected  to  become  available  shortly. 

. 

(2)  Study  of  the  EMI  Internal  to  IC  Op  Amps 

The  program  NCAP  calculates  the  nonlinear  transfer  function 
Hn(f|»...tf  )  by  first  calculating  the  circuit  admittance  matrix 
Y(fj+...+fn)  and  the  n-th  order  source  vector  In(f^,...,fn)^61^ 
and  solving  the  equation 


[Y(f 


^+. .  .+fn)]Hn(fi « . . .  »fjj) 


Vf 


i . V 


(5-1) 


It  may  be  possible  to  study  the  behavior  of  EMI  generation  inside 
an  Integrated  circuit  by  examining  the  admittance  matrix  and 


current  source  vector.1  J  An  alternate  method  was  suggested  in 
Section  4.2.  (See  p.  85.) 

(3)  Op  Amp  Modeling  Below  1  MHz 

All  efforts  to  predict  EMI  in  IC  Op  Amps  have  not  been  very  ac¬ 
curate  at  frequencies  below  1  MHz.  This  includes  predictions 
for  the  second  order  nonlinear  transfer  functions  using  the  741 
full  model  and  macromodel  and  the  LM10  macromodel.  Part  of  the 
problem  Is  a  result  of  the  very  low  EMI  observed  at  RF  frequencies 
below  1  MHz.  While  EMI  at  low  RF  frequencies  may  be  of  less 
practical  concern,  research  Into  low  frequency  modeling  of  EMI 
In  Op  Amps  may  contribute  to  a  better  understanding  of  Op  Amp 
full  models  and  macromodels. 

(4)  Macromodels  For  Linear  Electronic  Systems. 

Electronic  systems  designed  to  be  linear  usually  contain  non¬ 
linear  devices  such  as  BJTs  and  FETs  in' which  EMI  effects  may 
occur.  These  systems  often  contain  a  cascade  of  linear  stages 
and  a  variety  of  feedback  paths.  It  would  appear  that  a  general 
nonlinear  macromodeling  procedure  could  be  developed  for  such 
systems.  The  simplest  procedure  would  be  one  based  upon  the 
approach  used  by  Boyle  et  al.  The  input  stages  could  be  modeled 
In  a  physical  manner;  interior  gain  stages  and  output  stages  could 
be  modeled  as  linear  systems.  A  more  systematic  approach  would 
be  to  use  the  cascading  relationships  given  in  Chapter  4,  especially 
when  Internal  feedback  Is  used,  to  estimate  what  stages  make  important 


contributions  to  different  nonlinear  transfer  functions.  Stages 
that  can  make  an  important  contribution  to  the  nonlinear  transfer 
function  of  Interest  should  be  represented  by  an  appropriate  small 
signal  nonlinear  model.  Other  stages  can  be  represented  by  small- 
signal  linear  models. 

* 

(5)  Nonlinear  Macromodels  For  Digital  IC's. 

Macromodeling  of  digital  ICs  preceeded  that  for  IC  Op  Amps  and  is 
widely  used  by  the  semiconductor  Industry.  Unfortunately,  the 
digital- IC  macromodeVs  developed  appear  to  be  proprietary.  This 
may  mean  that  EMC  analysts  will  have  to  develop  their  own  digital 
IC  macromodels.  It  is  anticipated  that  small-signal  nonlinear 
macromodels  and  the  computer  program  NCAP  cannot  be  used  to  pre¬ 
dict  EMI  in  digital  ICs.  Large-signal  macromodels  and  time  domain 
analysis  routines  such  as  the  transient  analysis  routine  of 
SPICE2[53^  will  be  required. 
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MACROMODEL  PREDICTIONS  FOR  EMI  IN  BIPOLAR  OPERATIONAL  AMPLIFIERS 
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SUWARY 


A  aacromodel  containing  two  transistors  can  accurately  predict  how  amp¬ 
litude  modulated  RF  signals  with  RF  carrier  frequencies  in  the  .05  to 
100  MHz  range  are  demodulated  in  a  bipolar  operational  amplifier  (uA741 
op  amp)  to  cause  undesired  low  frequency  responses  related  to  the  modula¬ 
tion  envelope  of  the  RF  signals.  The  original  op  amp  macromodel  which 
was  developed  by  Boyle  et  al.  was  modified  by  adding  four  capacitors 
'(each  4  pF)  to' account  for  parasitic  capacitance  effects  in  the  integrated 
circuit  op  amp.  The  macromodel  can  be  used  when  the  RFI  signals  are 
Incident  upon  the  op  amp  input  terminals.  The  use  of  the  macromodel 
leads  to  a  substantial  reduction  in  computer  time  and  expense. 


1.  Introduction 

The  Nonlinear  Circuit  Analysis  Program  NCAP  provides  the  EMC  community 
with  a  usable  procedure  for  analysis  of  electronic  circuits  which  must 
operate  in  an  EMI  environment  (ref. 1-6).  Among  the  nonlinear  effects 
that  NCAP  can  predict  is  demodulation.  The  use  of  NCAP  to  predict  de¬ 
modulation  RFI  effects  in  linear  bipolar  integrated  circuits  (ICs)  such 
as  broadband  cascode  amplifiers  and  operational  amplifiers  (op  amps) 
has  been  demonstrated  recently  (ref. 7-8).  Up  to  now  all  NCAP  simula¬ 
tions  reported  upon  have  used  an  approach  in  which  every  transistor  is 
modeled.  For  example  the  25  transistors  shown  in  the  op  amp  circuit 
In  Fig.  1  were  each  modeled  in  ref .8.  The  objective  of  this  paper 
Is  to  demonstrate  that  the  op  agip  macromodel  shown  in  Fig.  2  which  con¬ 
tains  just  two  transistors  can  predict  accurately  how  amplitude  modulat¬ 
ed  (AM)  RF  signals  with  RF  carrier  frequencies  in  the  .05  to  100  MHz 
range  are  demodulated  in  741  op  aiftps  to  cause  undesired  low  frequency 
responses.  An  important  point  is  that  the  computer  time  required  for 
op  amp  macromodel  calculations  is  approximately  one- tenth  the  time 
required  for  op  amp  complete  model  calculations. 

This  paper  Is  organized  in  the  following  manner.  In  Section  2  the  op 
amp  macromodel  used  with  NCAP  is  shown.  In  Section  3  the  unity  voltage 
gain  op  amp  circuit  (voltage  follower)- used  in  the  investigation  is 
described.  Also  described  are  equivalent  circuits  using  the  op  amp 
complete  model  in  which  25  transistors  are  modeled  and  the  op  amp 
macromodel  in  which  2  transistors  are  modeled.  NCAP  predicted  RFI 
effects  are  compared  to  experimental  RFI  effects  in  Section  4.  The 
comparison  will  show  that  demodulation  RFI  effects  predicted  using  the 
op  amp  macromodel  are  in  good  agreement  both  with  those  predicted  using 
the  op  amp  complete  model  and  with  those  determined  experimentally. 
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Fig.  1  Schematic  of  yA741  operational  amplifier  (op  amp). 


2.  Op  Amp  Models 

The  computer  program  NCAP  contains  models  for  passive  components  such  as 
resistors,  capacitors,  and  inductors  and  for  semiconductor  devices  such 
as  pn-junction  diodes,  bipolar  junction  transistors,  (BJTs)  and  field- 
effect  transistors.  To  use  NCAP,  an  equivalent  circuit  for  the  electronic 
circuit  to  be  analysed  must  be  developed  using  the  available. models.  Then 
values  for  the  model  parameters  must  be  specified.  A  schematic  of  the 
uA741  op  amp  used  in  this  investigation  Is  shown  in  Fig.  1.  The  UA741 
op  amp  contains  25  BJTs.  The  direct  approach  is  to  model  the  op  amp 
completely.  Each  of  the  25  8JTs  in  the  op  amp  is  replaced  by  an  appro¬ 
priate  BJT  model.  Predictions  obtained  using  the  direct  approach  on 
how  amplitude  modulated 'RF  signals  are  demodulated  in  op  amps  to  cause 
undesired  low  frequency  responses  have  been  reported  upon  recently  (ref. 
8).  When  every  transistor  in  the  op  amp  is  modeled,  computer  simulations 
of  op  amp  circuit  performance  can  require  long  computation  times  and  can 
be  expensive.  An  alternative  approach  is  to  use  the  Integrated  circuit 
op  amp  macromodel  developed  by  Boyle  et  al.  which  can  accurately  predict 
op  amp  input  voltage  -  output  voltage  relationships  for  large-signal  non¬ 
linear  transient  operation  (ref. 9). 

Shown  In  Fig  2  Is  a  circuit  diagram  for  the  op  amp  macromodel  reported 
upon  by  Boyle  et  al.  Since  the  macromodel  is  to  be  used  with  the  com¬ 
puter  program  NCAP  which  performs  nonlinear  incremental  analyses,  sev¬ 
eral  modifications  have  been  made.  The  power  supply  terminals  have  been 
grounded.  Four  diodes  used  to  provide  current  limiting  and  voltage 
saturation  have  been  omitted  because  they  are  normally  reverse-biased. 

The  last  modification  involves  adding  capacitors  CJub  from  collector  to 

ground  and  from  emitter  to  ground  for  both  macromodel  transistors.  The 
capacitors  CJub  are  added  to  account  for  parasitic  capacitances  associat¬ 
ed  with  the  pn-junctions  used  to  provide  Isolation  between  transistors 
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Fig.  2  Circuit  diagram  of  op  amp  macromodel  as  notified  for 
us«  with  Nonlinear  Circuit  Analysis  Program  NCAP. 

Not*  th*  Inclusion  of  four  parasitic  substrate  capac- 

<UnCM  Csub- 


In  Integrated  circuit  op  amps.  The  importance  of  accounting  accurately 
for  parasitic  junction  capacitances  in  Integrated  circuit  op  amps  was 
discussed  previously  when* the  op  amp  complete  model  was  used  to  make  RFI 
predictions  (ref .8).  It  will  be  demonstrated  in  Section  4  that  It  is 
•Iso  Important  to  account  for  the  effects  of  parasitic  Junction  capaci¬ 
tances  when  the  op  amp  macromodel  Is  used  to  make  RFI  predictions. 

The  op  amp  macromodel  shown  In  Fig.  2  contains  two  BJTs.  The  NCAP  model 
for  th*  BJT  Is  the  nonlinear  T  model  (ref. 1-6).  The  nonlinear  T  model 
parameter  values  for  the  two  BJTs  in  the  op  amp  macromodel  are  given  in 
Table  1.  It  Is  appropriate  to  comment  briefly  upon  these  values.  The 
two  transistors  Q1  and  Q2  In  the  op  amp  macromodel  are  in  a  differential 
pair  configuration  that  corresponds  to  the  differential  pair  formed  by 
Q1  and  Q2  in  the  op  amp  schematic  shown  in  Fig.  1.  The  NCAP  model  para¬ 
meter  values  for  all  25  BJTs  in  the  op  amp  are  given  in  Table  1  of  ref. 8. 
The  NCAP  parameter  values  for  transistors  Q1  and  Q2  in  the  op  amp  com¬ 
plete  model  are  also  used  for  transistors  Q1  and  Q2  in  the  op  amp  macro- 
model.  These  are  two  exceptions:  (1)  a  dc  collector  current  Ig  >  10  uA 

for  Q1  and  02  in  the  op  amp  macromodel  Is  used  because  that  value  is 
specified  by  Boyle  *t  a 1 . ;  (2)  the  values  for  the  parameter  a  were  ad¬ 
justed  to  obtain  the  values  for  th*  dc  common-emitter  short-circuit 
current-gain  h„.  and  h__,  specified  by  Boyle  ct  al .  (ref. 9). 


TAKE  1 


NCAA  INPUT  PARAMETERS  FOR  uA7«l 
OP  AMP  NACROMOOEl 


Parameter 

ni 

02 

« 

4.34 

4.34 

*c  <*) 

14.S7 

I4.S7 

20 

20 

» 

0.1CS 

o.iss 

Iq  w 

10.C-4 

10.E-4 

O.SC-1 

Q.SE-3 

a 

0.90 

0.(2 

Vtma* 

400 

400 

k(F-V'*J 

1.23E-I2 

1.23E-I2 

Ref. 

1.091 

1.091 

CU  (F) 

1.23E-12 

1.23E-12 

4  (F/A) 

9.09C-9 

9.09E-9 

rfc  (earn*) 

930 

<30 

re  (otam) 

S.33ES 

5.J3E6 

C,  (V) 

0. IE-12 

0. IE-12 

C,  <F) 

0.  IE-1 2 

0.1E-12 

The  situation  of  Interest  Involves  electromagnetic  radiation  (EUR) 
Incident  upon  a  system  outer  enclosure  being  coupled  through  apertures 
In  the  skin  to  the  system  interior.  The  interior  EM  fields  induce  RF 
voltages  on  the  system  cables.  The  RF  voltages  are,  conducteo  to  semi¬ 
conductor  devices  such  as  Integrated  circuits  (ICs)'  located  inside 
electronic  equipment.  The  RF  voltages  can  cause  RFI  effects  in  ICs. 

The  specific  case  discussed,  in  this  paper  is  illustrated  in  Fig.  3. 

The  RF  voltage  induced  on  a  system  cable  Is  represented  by  an  RF 
signal  source  with  a  SO  a  Impedance.  The  bipolar  IC  is  an  operational 
amplifier  (op  amp)  which  Is  being  used  in  a  unity  voltage  gain  circuit 
called  a  voltage  follower.  The  RF  signal  which  is  an  amplitude  modu¬ 
lated  (AM)  sinusoidal  signal  is  being  injected  into  the  noninverting 
Input  terminal.  Previous  investigations  have  indicated  that  this  case 
exhibits  EM  susceptibility  at  lower  RF  power  levels  than  other  cases 
(ref. 10).  The  RFI  effect  investigated  in  this  paper  is  caused  by 
amplitude  modulated  RF  signals  being  demodulated  in  the  op  amp  to 
produce  undeslred  low  frequency  responses  related  to  the  modulation 
envelope  on  the  RF  signals. 

The  computer  program  NCAP  uses  an  iterative  procedure  to  calculate  non¬ 
linear  responses.  The  first-order  (linear)  responses  are  calculated 
first.  Then  the  first-order  node  voltages  are  used  to  calculate  second- 
order  current  generators  which  may  be  viewed  as  the  cause  of  the  de¬ 
modulation  being  investigated.  In  general  discrepancies  between  pre¬ 
dicted  and  measured  first-order  node  voltages  lead  to  larger  discrep- 


107 


Fig.  3  Schematic  Illustrating  the  unity  gain  op  amp  voltage 
follower  circuit  investigated.  The  desired  signals 
(unmodulated  CU)  and  the  EMI  signals  (AM  modulated 
RF)  were  connected  directly  to  the  noninverting  input 
terminal . 


ancles  between  predicted  and  measured  second-order  node  voltages.  For 
this  reason  considerable  effort  was  expended  In  developing  an  accurate 
linear  equivalent  circuit  for  tne  voltage  follower  circuit  shown  in 
Fig.  3  (ref. 8).  Shown  In  Fig.  4  Is  the  NCAP  coding  circuit  diagram 
for  the  voltage  follower  circuit  with  the  op  amp  complete  model  from 
ref .8.  Shown  in  Fig.  5  is  the  same  circuit  with  the  op  amp  macromodel 
'Substituted  for  the  op  amp  complete  model.  NCAP  predictions  made  for 
the  circuits  shown  in  Figs.  4  and  5  will  be  compared  to  measured  results 
in  the  next  section. 


4.  Comparison  of  Experimental  and  Predicted  RFI  Results 

The  purpose  of  the  investigation  is  to  determine  how  well  Che  op  amp 
macromodel  can  predict  the  demodulation  of  RF  signals  in  an  op  amp.  This 
will  be  accomplished  by  comparing  NCAP- predictions  made  using  the  op  amp 
macromodel  to  NCAP  predictions  made  using  the  op  amp  complete  model  and 
to  experimental  results.  The  experimental  results  were  obtained  using 
the  measurment  system  shown  in  Fig.  6.  The  op  amp  voltage  follower  was 
excited  by  amplitude  modulated  (AM)  sinusoidal  RF  signals.  The  RF  carrier 
frequency  fRF  was  varied  over  the  range  .050  to  100  MHz.  The  AM  modula¬ 
tion  Index  m  was  0.5.  The  AM  modulation  frequency  f ^  was  400  Hz.  The 
voltage  follower  output  voltage  at  400  Hz  was  measured  on  a  tuned 
audio  frequency  (AF)  voltmeter.  The  AF  voltage  VM  is  related  to  the 
voltage  follower  second-order  transfer  function  Hjlf^.-fj)  by 


Fig.  4  NCAP  coding  circuit  diagram  for  the  voltage  follower 
circuit  with  op  amp  complete  model.  The  circuit 
shown  corresponds  to  the  experimental  system  shown 
In  Flq.  6. 


Fig.  5  NCAP  coding  circuit  diagram  for  the  voltage  follower 
circuit  with  op  app  macromodel .  The  circuit  shown 
corresponds  to  the  experimental  system  shown  in  Fig.  6. 


*1  “W  W"  RF  frequency 

*2*W-fAF  f/^AF  frequency 

Modulation1  50%  AM  at 


Fig.  6  Experimental  system  for  measuring  the  second-order 
transfer  function  ^(f^.-fj). 


VM(rms)  -  (0.707)mA2H2(f1,-f2)  (1) 

where  m  Is  the  modulation  Index,  A  the  voltage  amplitude  of  the  RF  signal 
generator,  fj  the  RF  carrier  frequency  fRF,  and  f^  the  lower  sideband 

frequency  (fRF  -  f Aj;).  In  obtaining  £q.  (1)  It  was  assumed  that  the 

upper  and  lower  sidebands  contributed  equally  (ref. 11).  The  voltage 
amplitude  A  can  be  calculated  using 

*  ■  «> 
where  P^jn  Is  the  available  power  from  an  RF  signal  generator  with  Inter¬ 
nal  Impedance  R  .  Using  Eqs.  (1)  and  (2),  experimental  values  for  the 
gen 

second-order  transfer  function  H-(f. ,-f,)  can  be  determined  In  a  manner 
that  now  will  be  described.  *  * 


Shown  In  Fig.  7  Is  a  plot  of  the  tuned  AF  voltmeter  rms  voltage  reading 

»„  versus  the  generate*-  available  RF  power  9 _  For  9mmm  <  +5  dBm  the 

n  gen  gen 

VM  values  Increase  10  dB  for  each  5  dB  Increase  In  P  ,  and  Equation  (1) 
h  gen 

Is  valid.  For  Pg#n  >  +5  dBm,  the  values  begin  to  Increase  more  rapidly, 

and  Equation  (1)  Is  no  longer  valid.  Also  shown  on  Fig.  7  is  a  dashed 
line  which  corresponds  to  the  values  predicted  by  NCAP  for  second- 

order  nonl Inearl  ties  in  the  op  amp  voltage  follower  circuit.  Higher- 
order  nonlinear  effects  which  may  be  Important  for  PgefJ  »  +$  dBm  have 

not  been  considered  In  this  Investigation.  The  second-order  transfer 
function  H^fj.-f,,)  may  be  interpreted  as  a  normalized  nonlinear 

(demodulation)  response  for  the  op  amp  voltage  follower  circuit  at 
sufficiently  low  RF  power  levels  (i.e.  Pg#n  <  +5  dBm).  Experimental 

values  for  ^(fj,^)  can  be  determined  at  any  point  on  the  straight 
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Fig.  7  Measured  values  for  the  rms  audio  frequency  (AF) 
output  voltage  VM  vs  RF  Input  power  Pg#n  for  the 

op  asp  voltage  follower  circuit. 


line  drawn  through  the  data  for  Pgaf)  <  ♦$  dBm.  Using  Eq.  (1)  and  Eq.  (2) 
w1th  R  -  50  ohms  and  the  VM  value  at  Pgan  -  -10  dta.  experimental 
values^for  H^f,,-^)  were  determined.  These  values  are  given  in  Table 
2  and  are  plotted  In  Fig.  8. 


Also  shown  in  Fig.  8  are  values  for  M2(fJt-f2)  predicted  by  NCAP  for 

the  voltage  follower  circuit.  One  set  of  predicted  values  were  obtained 
using  the  op  amp  complete  model  shown  In  Fig.  4;  all  parasitic  substrate 
capacitances  were  assigned  values  of  2  pF.  The  other  set  of  predicted 
values  were  obtained  using  the  op  amp  macromodel  shown  in  Fig.  5;  all 
parasitic  substrate  capacitances  were  assigned  values  of  4  pF.  The  NCAP 
Input  coding  list  for  Fig.  4  Is  given  In. ref .8.  The  NCAP  input  coding 
list  for  Fig.  5  is  given  In  Table  3.  It  is  observed  that  the 

values  predicted  using  the  op  amp  macromodel  agree  well  with  the 
H2(fi,-f2J  values  predicted  using  the  op  amp  complete  model.  For  values 

of  H2(flt-f2)  >  -40  dB,  the  differences  are  less  than  5  dB.  For  values 
of  Hjtfj.-fj)  <  -40  dB  (a  low  value),  the  differences  increase  but  are 

still  less  than  10  d8.  The  agreement  between  the  op  amp  macromodel 
predictions  and  the  op  amp  complete  model  predictions  is  very  good.  The 
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TABLE  2 

SECONO-OROER  TRANSFER  FUNCTION  MEASUREMENT 
11*7*1  UNITY  GAIN  SUFFER  WITH 
Rgwi  •  -  10  dSa  *nd  •  *00  Hz 


A f  Freq. 

Af  vh  Read. 

WV* 

V<H*J 

VM  (rat  *V) 

(dB) 

SO  k 

0.008 

-44.94 

SO  k 

0.011 

-42.17 

70  k 

0.01S 

•59.48 

SO  k 

0.019 

-57.42 

SO  k 

0.023S 

-55.58 

100  k 

0.029 

•53.75 

ISO  k 

O.OS4 

-4  6.08 

200  k 

0.112 

-42 

2S0  k 

a. is 

-37.9 

300  k 

0.27 

-34.37 

400  k 

0.92S 

-28. 6 

S00  k 

0.92 

-23.724 

800  k 

l.S 

•19.48 

700  k 

2.2 

-16.15 

SOO  k 

3.0 

-13.44 

900  k 

3.8S 

-11.29 

1  H 

4.7 

-9.54 

i.s  n 

7.1 

-5.97 

2  H 

S.4 

-4.88 

3  H 

S.4 

•4.88 

4  N 

7.0 

-4.2 

S  H 

7.9 

-4.2 

S  N 

7.4 

-S.S 

7  N 

1.2 

-4.72 

B  H 

9.0 

-3.92 

9  H 

9.8 

-3.18 

10  N 

10. S 

-2.58 

IS  H 

13. S 

-0.393 

20  N 

13.S 

-0.393 

30  N 

14. S 

-0.23 

*0  N 

13.5 

•0.393 

SO  N 

12.0 

-1.42 

so  n 

11.0 

-2.17 

•  The  frequency  f,  •  1u, 

end  the  frequency 

•  fff  xNere 

fu  •  400  HZ. 

agreement  between  measured  and  predicted  results  shown  In  Fig.  6  may 
also  be  characterized  as  being  quite  good. 
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Shown  In  Fig.  9  are  three  sets  of  predicted  values  for  H^(f ^.-f^)  obtained 

using  the  op  amp  macromodel  shown  In  Fig.  5.  The  values  assigned  to  the 
parasitic  substrate  capacitance  C$ub  are  0  pF,  2  pF,  and  4  pF.  It  is 

apparent  that  the  value  assigned  to  Cj(jb  has  a  major  effect  upon  the 
Hjtfi.-fj)  values  predicted  as  shown  in  Fig.  9.  The  value  CJub  *  4  pF 
yields  the  best  agreement  between  predicted  and  measured  results. 


f, ,  RF  Frequency  (Hr) 


Fig.  8  Measured  and  predicted  values  for  the  second-order 
transfer  function  ^(f^-fj)  »$  RF  frequency  for 

the  op  aap  voltage  follower.  The  measured  values 
are  from  Table  2.  One  set  of  predicted  values  were 
obtained  using  the  op  amp  complete  model  shown  in 
Flo.  4  with  all  parasitic  substrate  capacitance 
values  set  at  2  pF.  The  other  set  of  predicted 
values  were  obtained  using  the  op  amp  macromodel 
shown  In  Fig.  5  with  all  parasitic  substrate  capac¬ 
itance  values  set  at  4  pF. 


5.  Conclusion 

A  macromodel  has  been  used  to  predict  RFI  effects  In  a  uA741  op  amp 
circuit.  Specifically,  the  computer  program  NCAP  was  used  to  calculate 
how  amplitude  modulated  RF  signals  incident  upon  the  Input  terminals  are 
demodulated  In  an  op  amp  voltage  follower  circuit  to  produce  undesired 
low  frequency  responses  related  to  the  modulation  envelope  of  the  RF 
carrier.  A  comparison  between  predicted  and  experimental  values  Indicates 
that  the  macromodel  can  be  used  quite  successfully  to  predict  the  de¬ 
modulation  RFI  effects  in  an  op  amp  circuit  for  RF  frequencies  up  to 
100  MHz.  The  op  amp  complete  model  can  be  used  for  RF  incident  upon  any 
op  amp  terminal.  The  use  of  the  op  amp  macromodel  is  restricted  to  the 
situation  In  which  the  RF  Is  Incident  upon  the  op  amp  input  terminals. 

This  usually  Is  the  most  susceptible  case.  When  the  RF  signals  are  In¬ 
cident  upon  the  op  amp  Input  terminals,  the  use  of  the  macromodel  (instead 
of  the  full  model)  leads  to  a  substantial  reduction  In  computer  time  and 
expense. 
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Curve  Csijb 
- 0  pF 

//  . **r 

•  //  -  4  PF 

—I - 1  I  I  .i.'l  f - 1 _ III  ■  .  .  ■ 

100k  IM  2  4  6  8 IOM 

f,,  RF  Frequency  (Hz) 

Measured  and  predicted  values  for  the  second-order 
transfer  function  H2(f]t-f2)  vs  Rf  frequency  for 

the  op  amp  voltage  follower.  The  measured  values 
are  from  Table  2.  The  predicted  values  were  obtained 
using  the  op  amp  macromodel  shown  in  Fig.  5  with 
all  parasitic  substrate  capacitance  values  set  at 
0,  2,  or  4  pF/ 


An  Important  point  is  that  the  op  amp  macromodel  developed  by  Boyle  et  al. 
has  to  be  modified  for  RFI  predictions  at  RF  frequencies  greater  than  1 
MHz.  Four  capacitors  CJub  are  added  to  account  for  parasitic  capacitances 

associated  with  the  pn-junctlons  used  to  provide  isolation  between  tran¬ 
sistors  In  the  Integrated  circuit  op  amp.  If  these  capacitors  are  omitted 
(CJ||b  *  0  pF),  the  RFI  effects  are  underestimated  by  more  than  10  dB. 

The  value  CJub  ■  4  pF  yields  good  agreement  between  predicted  and  measured 

RFI  results. 
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NCAP  COOING  FOR  VOLTAGE  FOLLOWER  WITH  HACROMOOEL 


•*  HACROMOOEL  of  741  OP  AMP 
•*  ANALYZED  GY  NCAP  FOR  H2 
**  SUBMITTED  AT  SONY  BUFFALO 
—  GORDH  OUTPUT  NOOE  IS  IS 
**  WITH  SUBSTRATE  CAPACITANCES 
“  ADJUSTED  PARAMETER  TO  SUIT 
•*  LM74I  HACROMOOEL  SETA  PARAMETER 

*  start  circuit 

•  GENERATOR 
NOOE  I  0 

FR  1  8.0ES  80.0E6  5  LIN 
AMP  1  0.0 

FR  2  -7.9996E5  -79.999ES  5  LIN 
AMP  1  0.0 
IMP  SO  0 

•LINEAR  COMPONENTS 
R  2  1  0.01 
R  2  3  560 
R  3  8  0.01 
R  2  0  SO 
*0  6  SMS 
R  7  LB  2712 
R  0  10  5305 


11  18  2712 
16  0  9.872E6 

12  0  l.OES 

13  0  42.87 

13  14  32.13 

14  IS  4.7E4 

15  0  4.489E4 
14  4  620 


C  2  0  1. IE-10 
C  6  10  S.460E-12 
C  16  0  2.41E-12 
C  12  13  3.0E-11 
C  IS  0  2.106E-9 
**  SUBSTRATE  CAPACITANCES 
C  6  0  4.0E-12 
C  10  0  4.0E-12 
C  7  0  4.0E-12 
C  11  0  4.0E-12 

•  TRANSISTOR 
NOOE  4 

4.34  14.57  20  0.165  10.QSE-6  0.SE-3 
0.90  400  1.23E-12  1.091  1.23E-12 
9.09E-9  830  S.33E6  0.1E-12  0.1E-12 

•  TRANSISTOR 
NOOE  8 

4.34  14.57  20  0.16S  10.0SE-6  0.SE-3 
0.62  400  1.23E-12  1.091  1 .23E-12 
9.09E-9  830  5.33E6  0.1E-12  0. IE-12 

•  LINEAR  DEPENDENT  SOURCE 
NOOE  C  10  0  12 

VC  1.886E-4  0.0 

•  LINEAR  DEPENDENT  SOURCE 
NOOE  16  0  12  0 

VC  6.23E-9  0.0 

•  LINEAR  DEPENDENT  SOURCE 
NOOE  12  0  0  13 

VC  247.49  0.0 

•  END  CIRCUIT 

•  ENO 
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Appendix  B 

IXLO  Op  Amp  and  Voltage  Reference  Data  Sheets 

National  0iC~ 

Semiconductor 


LM10/LM10B(L)/LM10C(L) 

Op  Amp  and  Voltage  Reference 


general  description 

The  LM10  seri«  an  monolithic  Ilnur  ICs  consisting  of  • 
precision  reference,  an  adjuitabia  rafaranc*  buffer  and 
an  indepvndent.  high  quality  op  amp. 

Tha  unit  can  opetsa  from  a  total  supply  voitaga  at  low 
at  1.1V  or  aa  high  at  40V.  drawing  only  270mA.  A  com¬ 
plementary  output  stage  savings  within  IS  mV  of  the 
supply  terminelt  or  anil  deliver  £20  mA  output  current 
wWs  ±0.4V  saturation.  Reference  output  can  be  as  low 
at  200  mV.  Some  other  cherecteriatiea  of  the  LM10  are 


•  bspueoffwt  voitaga 
■  inptit-offiet  current 


■  reference  regulation 

•  offsat -voitaga  dnft 

•  raParanea  drift 


10  mV  (max) 
0l7  nA  (mwi) 
a  nA  (max) 
0l  1%  (max) 
2|iV/*C 

0.002VC 


Tha  circuit  it  recommendad  fer  ponabta  aquipmant  and 
it  eomptataly  soaofiad  for  optration  from  a  fingla  power 
Mil.  In  contraat  hipt  output-drive  capability,  both 
voitaga  and  current,  dong  with  thermal  overioad  pro- 
action,  auggaat  it  in  damanding  ganer at -purpose 
applications. 

The  dance  it  capable  of  operating  in  a  floating  mode, 
independent  of  fixed  supplies.  It  can  function  as  a  remote 
comparator,  signal  conditioner.  SCR  controller  cr  trans¬ 
mitter  for  analog  signals,  dtlivaring  the  processed  signal 
on  tha  same  line  used  to  supply  power.  It  it  also  witad 
for  operation  in  a  wide  range  of  voitaga-  and  currant- 
regulator  aooiications.  from  low  vdtagat  to  savaral 
hundred  volta,  providing  gr attar  precisian  then  existing 

iCs. 

This  series  it  aeailaWe  in  tha  thru  standard  ttmparaturt 
ran  gat,  with  the  commercial  pan  having  r  at  axed  limits. 
In  addition,  a  low-vottaga  specification  (suffix  "L")  it 
aveilibis  in  me  limited  ttmparatura  ranges  at  a  coat 


d  teas  rwnv  tiwuniarnn  Cere. 


iMoafMiaa/vrwtve  m  us.  a. 


.  -  V-V- 

-v-T -U~  -- t-  .  V*-. ~  *.•  v  vl  s. 


LM10/LM10B(L)/LM10C(L)  Op  Amp  and  Voltage  Reference 


•lectrical  characteristics  it  j  -  2s*  c,  t**  £  t4  <  t***  not*  4) 


(liMfM  typt  iifwi  tm  Umio  ow  tampmom  raw  pi.) 


mmamctm 

CONDITIONS 

1  LM108L 

i  LM10CL 

UNITS 

MIN 

TYP 

MAX 

MIN 

TYP 

MAX 

Input  ottaat  sol  tap" 

13 

10 

0.S 

4.0 

mV 

3j0 

S.0 

mV 

Input  oOfaateunant 

ai 

0.7 

12 

2.0 

wA 

Inoat  SI 

IS 

10 

nA 

Input  Moo  ourtottt 

HI 

20 

12 

30 

nA 

30 

40 

nA 

Iftpyt  rwranci 

250 

500 

ISO 

400 

kfl 

ISO 

IIS 

kQ 

Ltrpa  sspnal  nlttpt 

V, -tl2SV.  lour  -0 

•0 

300 

40 

300 

V/mV 

VoefT  '  t!2V 

40 

20 

V/mV 

V«  •  H2SV.  Ipur  •  10  mA 

to 

2S 

s 

28 

V/mV 

VQUTa*2-7*V 

4 

3 

V/mV 

V.-tOdV  I0.S8V1.  loyr-tJmA 

1.S 

10 

1.0 

10 

V/mV 

Vour  •  taov  (tojvi.  Veto  •  hmv 

0.5 

0.7S 

V/mV 

Stunt  ftm  Inoaa  4) 

i.sv  <  v*  <  «.sv,  •  soon 

s 

30 

6 

30 

V/mV 

- 

0.1  mA  <  Iqut  <  10  mA 

4 

4 

V/mV 

Common  moda 

-3.25V  <  Vqu  <  2.4V  I2JSV) 

so 

103 

so 

102 

dO 

rajaedd* 

VftlSSV 

S3 

KB 

Shpphftnltut 

-02V>V*>-5.«V 

SO 

SO 

so 

90 

(» 

rawed  ow 

V*-IJV  C1.2V1 

SO 

d0 

1.0V  (1.1V)  ^  V*  ^  SJV 

04 

too 

so 

100 

dfl 

V-ajv 

m 

HO 

(Mat  rutupa  drift 

- 

10 

so 

uv/*c 

Offaat  currant  drift 

10 

S.0 

pA/*C 

•  Plot  ourrant  tpift 

« 

so 

pA/*C 

Una  topdation 

1-2V»1JVI<VtS«.SV 

02101 

0.01 

0.001 

0.02 

%/V 

OS  'npnS0.SmA.Vndn  -  200  mV 

0212 

0.03 

VV 

Load  mpdaddn 

OS'nsnS  O-StnA 

0.01 

0.1 

0.01 

ns 

% 

V  -  Vftdn  >  I-0V  (l.ivi 

IIP 

12 

% 

AmpNflar  pmn 

OJVSVn tnSUV 

30 

70 

20 

70 

V/mV 

10 

IP 

V/mV 

Pppdbaak  «aw«a 

too 

200 

206 

ISO 

200 

210 

mV 

ntum 

IPS 

200 

ISO 

211 

mV 

NaM  currant 

c. 

20  ' 

SO 

22 

70 

nA 

SP 

SO 

nA 

MaParanac  drift 

0.002 

0.003 

VC 

SupMv  currant 

200 

400 

200 

500 

nA 

soo 

570 

nA 

Nm  It  71w  input  wNtafi  can  «mm4  aw  wppiy  w*wpw  prpyipcd  dwt  aw  volte*  ,r<*"  aw  wove  to  toy  ottwr  tarmwwi  don  not  taend  aw 
nww wiww  awtramg  iww  Wtn  and  cuons  aiwpaaow  •»  tccognwd  far  ww  vt*  <  V~. 

Maw  2s  Tlw  wi— www,  oowatiwp  junction  awwwtwt  *  tSO*C  for  aw  UM10.  IOO*C  for  tfw  LMIOdIU  end  W*C  for  tfw  LM10CIL).  At 

otavawd  mmporaoumo,  dovwos  wwtt  bo  dorawd  bawd  ow  pecko*  dwrmai  rouwowco. 

Maw  2s  fwwrwd  awrwwa  i  muting  prevents  wctiwi  hooting  mat  could  result  w  sudden  faduro.  Out  dw  iC  can  ba  tubwetad  to  eceetoreted  *  tress 
wtdi  adwrwd  output  «wd  -ant  caw  cowdHWwo. 

NM*  T)m  ninmum  aoMv  tor  V-  <  Vcm  <  V*  -  OBSV  (lOVI,  11V  I1JVI  <  Wj  <  VMA*.  Vnfn  •  0.2V  w«0<  Iftfn  11  0mA. 

MOM  MMmM  MlfiM.  Vmax  *  <0V  far  ow  itandard  oan  and  6.5V  for  tfM  low  roKo«o  port.  Normal  tvoafaca  indicant  25* C  iimm.  PaMtata 

two  MMM  Mm  and  afaarao  Mat  tandltnnt  for  fiid-nmairanitt  ranoa  aoaranan:  Km  it  -55' C  to  >25*C  for  tot  CM  10.  -2S  C  to  t$*C  'or  tna 

UtlOWU  an*  ore  to  70*C  lor  tfw  LM10CICI.  To#  toaaficatiom  do  not  inctuda  ma  affaott  of  KMrmal  araotatm  Irt  l  20  mu.  «m  Hoatmo 

<T]  •  020  or  oatinoa  ftaotwQ.  OttdMnt  affaets  ara  mu  and  land  to  offtat  tfw  titctncai  arror  Itaa  cunm  i 

MaM  it  Oor  Tj  >»crc,  iQd  may  t teat <  UnA  lor  Vpu  -  V  nntn  Tj  -  125*e  nsd  V"  <  VCM  <  V- ♦  0.1V,  IQS  <  5  nA. 

data  0t  TMa  OaAnaa  aoaranan  in  ftootmo  aoonaanont  men  aa  KM  ooatttraepaa  rayuMtor  or  twomrira  trantmnttr.  Outout  >■  aonnactad  to  km 

V*  wrmlnal  of  KM  IC  and  <nowt  common  moda  is  rafanod  to  V  I  tar  tvowai  sooi  stations)  llfact  of  lar oar  outaumoitaoa  tmmoa  sntn  itifnar 

loodftdfnotctnut  ttadunttd  forut  tddlm  tnaraaittna  tuaoly  inaction  arror. 


mmrmmumm  mmimnmm 


VMMtVIfcptt 


typical  performance  characteristics  (op  amp) 


definition  of  terms 


Input  o that  wlmi!  That  voitagt  which  mutt  be  moiled 
between  the  input  term  inala  to  Out  the  unioadM  output 
in  tho  linoor  region. 

Input  oifnt  currant:  The  difference  in  tho  currant) 
at  tho  input  tasminaia  whon  tho  output  it  unloaded 
in  die  linear  region. 

Input  Hot  aunune  Tho  tbtohito  value  of  tho  at rarap 
of  dM  two  input  currant). 

Input  raUrtuiif  Tho  ratio  of  tho  ehtngo  in  input 
Niup  to  tho  chon  go  in  input  currant  on  ttthor  input 
with  the  other  grounded. 

Langs  stpul  vottego  gem  :  Tho  ratio  of  tho  specified 
output  voltage  swing  to  tho  ditngo  in  difftrantial 
input  voltage  required  to  produce  it. 

Shunt  goin:  The  ratio  of  tho  spoeifiod  output  voitagt 
swing  to  tho  change  in  differential  input  voitagt  required 
to  produce  it  with  tho  output  tied  to  tho  V*  terminal 
of  tho  1C  The  load  and  power  source  are  connected 
between  tho  V*  and  V~  terminals,  and  input  common* 
mode  it  referred  to  tho  V-  terminal. 

Common  modo  refaction:  Tho  ratio  of  tho  input  voitagt 
range  to  tho  change  in  offiot  voltage  Between  tho 
extremal. 


Supply  voltage  roioettom  The  ratio  of  tho  tpodfied 
supply-voltage  change  to  tho  change  in  offset  voltage 
between  dto  extremes. 


Lino  regulation:  The  average  change  in  reference  output 
voltage  over  the  specified  suppi\«  voltage  range. 


Lead  regulation:  The  change  in  reference  output  voitagt 
from  no  load  to  diet  load  specified. 


ftedberk  settee  voltage:  The  voltage,  referred  to  V” 
on  tho  reference  feedback  terminal  while  operating 
in  regulation. 


Wife  ranee  emptlfler  gain:  The  ratio  of  tho  specified 
reference  output  change  to  the  change  in  feedback 
Sanaa  voitagt  required  to  product  it 

Feedback  curtert  The  absolute  value  of  the  current 
at  the  feedback  terminal  when  operating  in  regulation. 

Supply  ourrant:  The  currant  required  from  the  power 
source  to  operate  the  amplifier  and  reference  with 
their  outputs  unloaded  and  operating  in  the  linear 
range. 
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|  MISSION  j 

§  ftww  i4tr  Development  Center  | 

ft  RAPC  p-tanA  and  execute*  ruearch,  de.veZopme.nt,  tut  and  r 

>  &  eluted  acquisition  programs  In  support  o^  Command,  Control  Q 
&  Communications  and  Intelligence.  (CSJ)  activitiu .  Technical  ? 
**  and  engineering  support  utithin  areas  of  technical  competence  C 
$  Is  provided  to  ESP  Program  Offlcu  l POs )  and  other  ESP  * 

u  elements.  The  principal  technical  mission  areas  are  % 

?  communications,  electromagnetic  guidance  and  control,  sur-  .  <* 
ft  veiUance  of  ground  and  aerospace  objects,  Intelligence  data  r 
\  collection  and  handling,  Information  system  technology,  Q 

fe  Ionospheric  propagation,  solid  state  sciences,  microwave 
%  physics  and  electronic  reliability,  maintainability  and 
fc  compatibility. 


